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EDITORIAL 


Beginning with this issue the Proceedings will appear in two Parts, 
“A” and “B”. ‘This change has been dictated by the increasing 
volume of important communications which has proved too great even 
for the larger amount of space available since monthly publication began 
in January 1948. 

It is intended provisionally to publish in Part A papers dealing 
with subjects such as the quantum theory, statistical mechanics, 
nuclear physics and cosmic rays, atomic physics, molecules, spectra, 
theories of solids, liquids and gases, surface phenomena, growth and 
properties of crystals, crystal structure, luminescence, electrodynamics, 
heat and thermodynamics, and standards. Part B will include subjects. 
such as acoustics (including ultrasonics), optical design, electron 
optics, colour, elasticity and other mechanical properties of solids and 
liquids, crystal structure analysis, magnetic materials, refrigeration and 
liquefaction, electric discharges, radio, geo- and ionosphere physics, 
astrophysics, and solar physics. It is obvious that a hard and fast 
dividing line cannot at present be drawn, but we feel sure that in a short 
time each Part will develop a character of its own, so that the reader 
will have little difficulty in deciding with which he wishes to be supplied. 


With this new arrangement we hope to be able, not only to give 
improved service to Fellows and other subscribers, but also to effect 
that speedy publication of important scientific papers which is so vital 


to the development of physics. 
G. INGLE FINCH, 


President. 


PROC. PHYS. SOC. LXIT, I—-=A 


Collision Problems and the Theory of Radiation Damping * 
By J. HAMILTON 


Manchester University 
Communicated by L. Rosenfeld; MS. received 12th May 1948 


ABSTRACT. The relation between radiation damping theory and the problem of the 
radiationless collision of an electron with a field of force is discussed. It is shown that 
the divergence difficulties, such as appear in the damping treatment of Compton scattering 
do not arise, and the ‘“‘ cut-off’ of damping technique is not necessary. The same result 
applies to the relativistic collision of two electrons. 


§1. INTRODUCTION AND SUMMARY 


T has been pointed out by Bethe and Oppenheimer (1946) that an uncritical 

application of the equations developed by Heitler and Peng (1942) for 

treating radiation damping, leads to incorrect results when applied to the 
radiationless collision of particles under a static interaction potential. It is, of 
course, clear that the damping theory was developed in order to treat the inter- 
action of quantized fields with particles to a higher degree of accuracy than is 
possible by the conventional perturbation theory. It may, however, be of interest 
to discuss the relation between damping theory and collision problems from a 
general point of view. This should give some idea of the range of validity of 
the theory, and some guide as to its application to more complicated problems, 
such as bremsstrahlung. 

The scattering of a particle by a central static field of force is considered 
in the non-relativistic case. ‘This collision problem is treated by the same method 
as was used in an earlier paper (Hamilton 1947), to be referred to as I, to deal 
with the damped scattering of radiation by an electron. ‘Thus it is possible 
to compare in detail the treatment of the collision problem and of Compton 
scattering. ‘The equations which describe the change in the free particle system 
due to the presence of the potential energy can be written in the same form as 
the equations which connect the initial and final states in the damping theory 
treatment of Compton scattering. The first difference between the two cases 
is that infinite transverse self-energy terms have to be neglected in the Compton 
scattering problem before any reasonable equations connecting the initial and 
final states can be formulated. 

A more important difference arises in the solution of these two sets of equations. 
In Compton scattering it is essential to limit the energy of the final states to a 
narrow band of values centred on the energy of the initial state. The excluded 
final states would give an infinite contribution of such a type that the scattering 
cross-section would vanish (cf. I, §7(a)). Relativistic invariance is not violated 
by this procedure, as the breadth of this energy band tends to zero as the size 
of the enclosing box becomes infinite. If the final states to be considered in the 
collision problem be limited to a similar narrow band of energies it is shown 
that a result is obtained which is in agreement with the result got by applying 


‘ “This paper forms part of a thesis presented to the University of Manchester for the Ph.D. 
egree. ; 
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the Born approximation, when the latter is valid. When the Born approximation 
is not valid, this limitation to a narrow band of final energy values is easily seen 
to be incorrect. 

Further investigation of the collision problem shows that those final states 
which lie outside the narrow energy band do not give an infinite contribution; 
in fact they give just that contribution which is necessary to lead to the correct 
collision cross-section as deduced by the phase method. Thus, in general, a 
direct application of the equations which were developed to treat the damped 
scattering of radiation leads to results which have no greater accuracy or range 
of validity than those given by the Born approximation. It seems that, when 
there is no reason (such as divergency) for the limitation of the range of energy 
values of the final states, the Heitler-Peng equations, in their simple form, are 
no longer valid. 

In §§ 2-4 the case of a central force of finite range is treated, and in §5 the 
modifications which are necessary in dealing with the Coulomb potential are 
indicated. In both these cases it is possible to simplify the calculations by using 
a simple relation which exists between the phase shift of the wave mechanical 
method of treating scattering, and the energy shifts which play an important 
part in the damping phenomenon. 

In §6 the relativistic collision of two electrons (without the emission of 
radiation) is discussed. While the restriction to radiationless collisions makes 
the discussion of limited validity, it seems that there is no more serious convergence 
trouble for high energies of the final states when the lowest order relativistic 
interaction is used, than when the Coulomb interaction alone is used. Thus the - 
parts of the interaction due to longitudinal and transverse photons, respectively, 
are not treated by different methods, as would be the case if the transverse part 
had to undergo a cut-off procedure; and the relativistic invariance is not violated. 

It would be useful to have a general criterion which would indicate those 
‘cases in which it is essential to limit the possible energies of the final states. It 
is a fact that, so far as the electromagnetic field is concerned, this limitation has 
only been required in problems in which a change in the free photon configuration 
is considered (i.e. emission or scattering of a photon). This may be a general 
tule. 

It might be expected from a correspondence principle argument that no damp- 
ing problem in the usual sense arises in the collision of two electrons. In classical 
or quantum theory, damping expresses the reaction of the emitted or scattered 
entity on the emitter or scatterer. In the former, damping takes the form of 
certain reactive forces, while in the latter it can be regarded as arising from a 
re-normalization of the wave function of the initial state. ‘The classical theory 
of electron collisions introduces no damping effect of this nature, and the same 
may be expected to hold for the non-relativistic quantum theory treatment. 

That divergence troubles are absent from the simple treatment of the radiation- 
less relativistic collision of two electrons is also to be expected. Moller (1931) 
has shown that the part of the interaction which results from the virtual emission 
and absorption of one photon, can be regarded as being due to the interaction of 
the currents which arise from the motion of the electrons. As it is possible to 
treat the motion of an assembly of charges and currents by the ordinary quantum 
mechanical methods without any difficulty, it is to be expected that no divergence 
of energy integrals will arise in the present treatment. It may be of interest 
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to remember that the relativistic interaction of the electrons can be regarded, 
by the correspondence principle, as arising from uniformly moving electrons, 
so that no acceleration of the electrons is involved, and no radiation reaction 
arises. 

These arguments apply to the electromagnetic field only. ‘The interaction 
of the charges and currents arising from the meson fields of nucleons cannot be 
treated by the ordinary quantum mechanical methods, and divergence difficulties 
are expected to occur in treating collisions involving these particles. 

There is, of course, in principle, a distinction between the existence of diver- 
gence difficulties and the presence of a damping effect. Even if quantum electro- 
dynamics were well behaved so that no divergent integrals arose, it would be 
possible to express the damping effects of emitted or scattered radiation. In 
that case the Heitler-Peng formalism would have to be altered so as to allow for 
those final states which lay outside the narrow energy band. In view of the above 
discussion it is possible to go further and state that such a modification would 
be particularly necessary in cases where the interactions are strong, and the 
damping is, therefore, really important. 

The question of the effect of the emission of radiation on the collision of an 
electron with a field of force is not considered in the present paper. Bethe and 
Oppenheimer (1946) applied the Heitler-Peng equations to the bremsstrahlung 
process in which only one photon is emitted, and they found that the infra-red 
catastrophe appeared in a very troublesome way. As all the final states interact 
with each other in any damping treatment it is essential that all the physically 
important states are considered. 

Ferretti and Peierls (1947) pointed out that this imposed a restriction on the 
validity of Bethe and Oppenheimer’s treatment, due to the possible emission 
of many photons. They showed that for the single photon treatment to be 
valid, it is essential to enclose the system in a reflecting or periodic box of a size 
such that (e?//ic) ln (E/Pyin) <1, where p,,;, 1s the lowest photon energy allowed 
by the size of the box, while £ is the initial kinetic energy of the electron. Ferretti 
and Peierls showed that, if the size of the box be such that this condition holds, 
there is a damping correction to the radiationless scattering cross-section, which 
depends logarithmically on the linear dimensions of the box. They suggested 
that this was a meaningless result. 

The bremsstrahlung problem has been partially investigated. It has been 
shown that, in the single photon case, the damping correction to the cross-section 
for radiationless collisions of the electron is just the total cross-section for collisions 
in which a photon is emitted. Although this correction does depend on-the 
size of the box in the way Ferretti and Peierls suggest, it clearly has a very reasonable 
interpretation. A solution of the bremsstrahlung process for the case in which 
the box may be infinite appears to be very complicated. 


S725 AMSUS IPISLANSS) Sisuiln sh 


In this and the two following sections the collision of a particle with a central 
field of force of finite range is considered. Spherical polar coordinates are the 
most convenient. In order to show the connection between the phase shift 
of the wave function of the scattered particle and the change in the energy levels 
which occurs in going from the free particle system to the system in which the 
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potential energy acts on the particle, it is essential to enclose the particle in a 
large spherical box. This box is assumed to be impenetrable; its centre coincides 
with the centre of force, and its radius is denoted by R. 

The motion of a free particle in such a box is given by the wave function 


DCCs PP) Ops © aur hsores (1) 


where J and Y are Bessel functions and spherical harmonics, respectively, and 
C is a normalizing constant; p is the radial momentum, and h is taken as being 
unity. The asymptotic form of such a wave function is ®~r~'!sin (pr —Iz/2). 
The boundary condition requires pR—J/7/2=nz, where n is an integer. Thus 
the p values are equally spaced by the amount 


Mpa eet es Ser Se (2) 


This corresponds to an energy spacing, in the non-relativistic case, of 
AE =(p/m)(7/R), where m is the mass of the particle. 

When the problem of the scattering of a free particle by the potential energy 
V(r) is solved by the normal methods, the result can be expressed by the asymptotic 
form of the wave function, viz. Y~r-!sin(pr—/7/2+68,). If this scattering 
takes place within an impenetrable box of radius R a spectrum of p values is 
selected, given by p’R—In/2+8,=m7r. These new values of the momentum 
are related to the momentum values of the free particle by the equation 


(p’ —p)R SO Viet” Cais 6) 


If A is the new energy value corresponding to p’, and F the original energy, then 
(2) and (3) give 
7A CN = 88 WER es) a an Ae (4) 


§3. THE BORN APPROXIMATION 


As the potential V(r) is independent of angles and spin the coupling of the 
free states (1) which V introduces does not mix the angular momentum quantum 
numbers /, m. In what follows the case 7=0 will be treated. The method of 
treating damping theory developed in I will be followed, in so far as it is applicable 
to the scattering problem. Denoting the states of the free particle by the quantum 
number p, the equations giving the behaviour of the system when V(r) is introduced 
are (cf. I, eqn. (4)) 

(Gig Ag+ iggy a0. Ow dene (5) 
Fy 


E, and A are the energies of the free and the coupled system respectively, the 
A, are the normal amplitude functions, and H,,,= a rdr®,(r)V(r)®,,(7). 
®,,(7r) is given by (1) with /=0. 

Following the treatment of Compton scattering (I, §6) new amplitude 
functions are introduced by the relation V,=(A—£,)A,. Substitution in 
(5) gives 
Vy = 2 Ayy Vy [(A- Ey): CY eee oars (6) 

Pp 
If A lies in the region between E,,, and the next highest energy value, the sum in 
(6) can be split up, conveniently, into a sum over the 2N values of £,, situated 
symmetrically about £,,,, and a sum over the remaining £,, values. 
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The result of this separation gives 


at ny Vp" 5 App Vy’ 
"p a A — Ey 3 “Ep En — Ey 
where [2N] denotes the region about £,,,, and (E,) the remainder. The details 
of the transition from (6) to (7) are given in the Appendix. 
Equation (7) gives 
Vie H Ver cot (AE \n/ME} + —| spy ap: (8) 
Pp pp" p ?" ie as yo) Bp = Eg. © ee teay > 
where the second sum in (7) has been replaced by an integral. (p’’) indicates 
that the principal value has to be used at p’ =p”, while H’,,,, = RA, 

According to the damping theory technique for radiation developed in I, 
the second term on the right hand side of (8) has to be neglected. In the radiation 
case this term is infinite. Neglecting this term in the present case will lead to 
the results obtained by indiscriminate application of the Heitler-Peng equations 
to collision problems. 

Thus (8) gives V,=H'p»-Vy(m/p") cot (A—E,,..)7/AE}. Putting p=p", 
it follows that 

tan {(A—E,)r/AE}=(m[p A yyy we eee (9) 
or using (4), 
tan Sy = —(Mip JH pope ree (9a) 


Allowing for the fact that the normalized functions ®,(r) of (1) have a radial 
part given by ®,,(r) =(zp/R)*r*J,(pr), (9 a) gives 


Sar2m (* 


tand)= — na \ VON PNPrdr ese (10) 


on re-introducing h. 
This agrees with the usual Born approximation for 59, viz. 


Sar2m1 


3)= dr | : Vir(prPrdr, seas (11) 


which is valid provided 6,<1. 

It can be seen readily that (10) gives incorrect results if the right hand side 
is large. Consider a potential V(r) which is extremely strong where r <a, and 
which vanishes where r>a. Then the true value of 8, is approximately —(a/p), 
but the right hand side of (10) is, in general, certainly not —tan(a/p). 

Thus the technique used for the damping of radiation when applied to a 
collision problem gives a value for the phase which agrees with that derived by 
using the Born approximation, in cases where the Born approximation is valid; 
but it gives an incorrect result in cases where the Born approximation is not valid. 
It follows that the region of validity of the Heitler-Peng equations for collision 
problems is the same as that of the Born approximation. 


§4. THE EXACT: SOUUION 


It remains to show that it is possible to take account of the second term on the 


right hand side of (8), and that this equation then leads to the correct value of 
the phase. 
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Equation (8) can be written 


PINE oen Ventas “| Vibes cee (Bia) 


ao are 


where d is an eigenvalue to be determined. H’ is independent of R, the radius 
of the box, and so is ®/ (7) = R'®,(r). Actually @3(r)=4/2 sin(pr)/r. To find 
the solutions of (8a) it is useful to study the Bie per yes of the homogeneous 
integral equation 


Woe ie oe Wy Bas (12) 


Equation (12) will have a set of real eigenvalues yu,, and a corresponding set of 
solutions Wf. The functions W% are orthogonal, and if the « form a discrete 


set, the Wz can be normalized so that | dp'WiW3,=8,5. ‘Then 
Ep SW W2 ligne ore be (13) 


where the summation extends over all values of «. However, the definition 


(ce) 


fA | BL UNA eet Uy (14) 


remembering the orthogonal properties of the M)(r), suggests that « denotes 
the continuous spectrum, given by 7, and that the Wj are related to the ®,(7). 
If Wj = *7®,(r) =(2/77)* sin(pr) then 

eae on late !} =a(7—1') 


a a r+yr’ 


lim [ WiW) dp=7 lim 


P+o J P+>o 


where 8(r—7’) is the Dirac 5-function. Similarly, the other orthogonal relation 
holds, viz. 


lim W; Widr=pap). 2 Wee (15) 


R>o 
Now (14) can be written in the form 


wre. 
Hy =| 2 Pe aig ae (16) 


which is analogous to (13). It follows that 1/u,=7V(r). It can easily be verified 
that the Wf, w,, thus defined are the solutions of equation (12). 
To solve (8 a) substitute V,, = i dr a, W; (i.e. the solution of (8 a) is expanded 


in terms of the solutions of (12)). 
Using (16) and the orthogonality relations (15) this gives 


oe dip N inte) teint ee we os (17) 
Wy Wy : : 
where (r, r')= WW, My [% File “{ ”) dp’ ie Tae -. Using Ee =p*/2m gives 


a2 apcosir'p.) sin (7p )-- r<7' 
[7, if ] SS (2m|7p i sin (r'p’’) cos (rp’’) ft. : 
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and also (r, r’)=(2/m)sin(p’r)sin(p’7’). Further, substituting «,=y4,, (17) 
leads to the equation for «, 


a, = 2{X sin (p’’r) —(m/p’’) cos ("7)}| dr'a,V(r') sin(p’’r’) 


—2(m|p’’) sin(p"’r) | j dr'a,V(r') cos(p"'7’). ee 


From (18) it follows that if «, is everywhere finite then %=0. By differentiation 
of (18) it is readily seen that «, satisfies the equation 


d2u,|dr2+{(p"2—2mV(r)}%p=0. nace (19) 


As is well known the solution of (19), together with the condition «)=0, leads 
to the asymptotic form for large 7 


bw din (6°74 89), eee (20) 


where A tends to a constant as roo. The asymptotic form of (15), on the other 


hand, is 
a,~2{X sin (p’’r) —(m/p’’) cos (p"7)}} dra VE) pT), sea (21) 
0 
Comparing (20) and (21), it follows that cot 69 = —Ap’’/m. 

Thus cot 8) = —cot {(A—E,,)7/AE}, in agreement with the exact relation (4). 
Hence the additional term in (8 a) gives a finite correction to the Born approxima- 
tion—a correction which gives the exact result as deduced by the phase method. 

Equation (21) yields the further relation for the phase 


sin dy = — 2¢mip")-*{ dr b(r) V(r) sin (pn), mage ae (22) 


where y(r) is the finite solution of (19) which has the asymptotic form 
w~A sin (pr +69). Equation (22) can be verified by the usual methods. 

Finally, it should be noted that the same method could be used for values 
of J other than zero, with a slightly more complicated analysis. 


5) DHE COU LONBS BibT D 


The investigations of the previous sections can be carried through for the 
Coulomb interaction, with some slight modifications. It is well known that 
in this case the phase contains an additional factor (—Inr). This corresponds 
to the fact that H’,,, depends logarithmically on R, the radius of the enclosing 
sphere. It is therefore necessary to keep the radius R finite, but large. The 
connection between the phase shift and the energy difference, established in 
§2, then holds (provided the In R term is included as being part of the phase). 
The p values now form a discrete set, and some integrations are to be replaced 
by summations. 

In the Appendix it is shown that the transition from (6) to (7) can also be made 
in the Coulomb case. The last term in equation (8) is to be replaced by the 


summation 7 ApS yn(H'yy)/(Ey Ey) where H' py = [ '(7)V(1)®'y (rr? dr. 
Pp 
Further, (15a) becomes Aph W/W?’ =8(r—r'). The modified form of (12), 
p 
W,=pAp=H',,,W,, has the solutions W%. Equation (17) can be derived, 


with the modifications that the upper limit of the integration over r’ is now R, 
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while [r, 7’] is given by a summation over the p’ values. Equations (18) and 
(19) follow easily, and (21) is replaced by 


tp =2{A sin(p’’R) —(m/p"’) cos(p’’R)} [ dr'«,V(r’) sin (p’’7’). 
0 


From this it again follows that cot y= —cot {(A —E,,,)7/AE}. 
It is clear, therefore, that for the Coulomb interaction the damping theory 
‘cut-off is no more justified than it is for a static interaction of finite range. 


§6. THE RETARDED INTERACTION OF TWO ELECTRONS 


The result of the preceding section leads to the question of what part damping 
plays in the collision of two electrons when the retardation and velocity dependent 
effects are considered. It might be guessed that as the damping theory was 
developed for quantized fields, it will modify that part of the interaction which 
_ arises from transverse photons. On the other hand, it is difficult to see how such 
a treatment could avoid destroying the Lorentz invariant form of the total inter- 
action (cf. Meller 1931). 

Any investigation of the nature of the corrections to be applied to the first 
order perturbation method of treating the relativistic scattering of electrons 
must be mainly formal. Of the two reasons for this, the first is that the correction 
-due to the emission of radiation is of the same order as any correction to Moller’s 
(1932) collision cross-section. ‘The other reason is that the matrix elements 
are very strongly dependent on the angle of scattering, if the latter is small; 
so that it is apparently impossible to apply perturbation theory to calculate the 
higher order corrections. 

It is easiest to use a linear momentum representation (and a cubic periodic 
ox). The fundamental equations (5) become 


(E; —A)A,y + 2A, evAct — 0, eo akemeee: (23) 
el 


where c,c’ label the distinct energy levels, while /, J’ denote the directions of 
the momenta and the spins of the two electrons in the initial and final states 
respectively. H,., is the matrix element connecting the states (c,/) and 
(c’, l'), and is composed of the Coulomb interaction matrix element, together 
with the compound matrix elements arising from the exchange of one transverse 
photon between the two electrons. 

Substituting (A—£,)A,=V,,, and changing the summation over /’ into 
an integration, equation (23) can be written as 


Fy, evP(Eel’) Vey ; 
Peake 


where A lies in the vicinity of E;. p(E,l)AE,dl is the number of final states 
in the range (AE,, dl); AE, is the interval between two discrete energy levels in 
the neighbourhood of E,,, while = cot {(A — E;)7/AE;}. (In the linear momentum 
representation there is no simple relation between A and the phase.) 

In the damping theory treatment of Compton scattering the second term on 
the right hand side of (24) is neglected, because the integration over the energy 
diverges. It can easily be seen that, in the Compton scattering case, the compound 
matrix element corresponding to’H,, »7 behaves like (E,)* as E,—> ©. 


Vy=ad| dl Hy erp(Ed' Vat | died Seether let, (24) 
(Ez) 
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Suppose that in the problem of the collision of two electrons the initial 
momenta are p,, p,; and the final momenta are p,, p;. Then the Coulomb. 
matrix element is 


4re*[|pa—Pal?], (A=c=1), sean (25) 
while the term in e? due to the transverse photon field is (Heitler 1944, p. 98) 


Dp rg (Ea FEB] (E, + E;) —(E, +E) -2k 
va [ ( Ra Pe (fe ee 


E,, E,; E,, E, are the initial and final energies of the electrons, while 
k=|p,—p,|=|p;—P»|- Conservation of energy is not assumed in calculating 
(25) or (26). 

For k much larger than p,, p,, and the rest mass of the electrons, (26) takes 
the form (27e2/k?)[a%a”) +1] and, for very large final energies, E,~2k. It is 
therefore to be expected that if the matrix element (25) gives no divergence 
difficulty in the energy integral occurring in equation (24), the same will apply 
to (26). As it is not permissible to neglect all but a narrow energy band in the 
case of the Coulomb interaction, it seems unjustified to apply the damping theory 
cut-off to the interaction given by (26). However, the general remarks on the 
essentially formal nature of any easel of (26) at high energies should 
be remembered. 

Finally, it may be noted that substituting the left hand side of (24) in the: 
integral on the right hand side gives 


Vu=za | ar) aa + bes dE, far Racer ens \ eB Va 


(26) 


+(a term in H.H.H.): 


Here | dE, is a convergent integral, but fav gives rise to an infinity when /’’ 


is parallel to J, or to J’ (in the centre of gravity system). Such a substitution 
method of approximating to (24) is therefore impossible, because of the strong 
angular dependence of the matrix elements. 
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APPENDIX 


The transition from equations (6) to (7) for finite range interactions can be- 
justified as follows. Consider first the (E,,,) region, which is the region remaining: 
when the N terms higher and lower than Ev are excluded. 

On (E,) the term 1/(A—E ue can be written 


1 1 
= >——— +(A-E,y) —— -(,+-) 
A =< E , E ia Ee - CE y» Ey” ce ie 


_ Ax Ep? za 1 
21 OE? \B,. Ey) a 
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The order of magnitude of the ratio of the second term on the right hand side 
to the first is |(A—E,,.)/(E,—E,,)|; and on (E,,”) this is less than r/N, where 
(A-E,,)/AE=r, and r= O(1). Choosing N= KR%9, where 0<e<} and K 
is a ponstant. this ratio can be made negligible. Thus aie second, and similarly 
the higher order terms on the right hand side of the series can be neglected 
compared with the first term. With this definition of N, it follows that 
NAE= (ap"/m)R— ») tends to zero as R+>co. Thus the breadth in energy of 
the region [2N] about £,,, which contains the 2N energy levels tends to zero as. 
R>o. 
For the region [2] it is convenient to use the expansion 


0 
Ay Vy = App Vy + (Ey — Ey) an ow Vr) + AR MOOR OS (A), 
p 
It is assumed that H,,,,,V,,, does not vary rapidly on [2N]. Further 
By — Ee 1 
2)? = —-2N+(A—-E,) & 
[2.N] A-— 1 ( Sf rao En’ 


If N, and so R, is sufficiently large, the last term gives 

{(A —Eyy,)n/AB} cot (A —E,-)n/AE}. 
It is therefore negligible compared with the term SON The first term on the 
right hand side of (A) gives the first term.in (7) whose value is 

FD yp V pA] RAE) cot (A —E,)7/AE}, 


pp” 


where H “pp” = RE, 


pp’ 


H’ 5 is independent of R provided 
R{\® alt)V(1)®,(r)r* dr, 


converges as R>oo. ‘This is true when V(r) is a potential of finite range. RAE’ 
is independent of R. The second term on the right hand side of (A) gives a 
term in (7) 
0 2N oO 
a= pr AE, \SesS = ea 3+8) 
N aE. (Hop V iy) R op, Ht vp V7) =O(R?**) 
This can be neglected. Similarly it can be shown that the higher order terms. 
in (A) do not contribute to (7). 
When I(r) is of the form 1/r it is essential to keep R finite, but it can be seen. 
that by choosing R sufficiently large, the errors involved in going from (6) to 
(7) can be made as small as is desired. 
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ABSTRACT. It is shown that a careful examination of the concepts involved in radiation 
damping theory makes it possible to treat the propagation of energy from one atom to 
another through the electromagnetic field. An exact and continuous solution of the 
propagation problem is given. 


§1. INTRODUCTION AND SUMMARY 


T has been suggested by Ferretti and Peierls (1947) that the quantum theory 

of radiation damping leads to incorrect results when it is applied to a problem 

which involves the propagation of radiation. They consider the resonance 
‘interaction of two atoms through the medium of the electromagnetic field; and 
they suggest that a consequence of the damping theory is that energy is propagated 
from one atom to the other at a velocity other than the velocity of light. They 
reach this conclusion after making use of the ideas of ‘‘ roundabout transitions” 
and “‘intermediate”’ states which are involved in the application of damping theory 
to the scattering of radiation (Heitler and Peng 1942). The present paper is an 
attempt to show that there need be no difficulty in the propagation problem, 
provided these ideas are examined carefully; and, in addition, it is an attempt 
to obtain an exact solution for this problem. 

It is assumed that there are two identical atoms A and B, which can interact 
by emitting and absorbing photons. For simplicity the states considered are: two, 
in which one atom is in its excited state and the other is in its normal state; and 
an infinity, in which both atoms are in their normal states and one photon is 
present. For such a system there need be no confusion due to ‘‘roundabout 
transitions”. This type of transition can be defined as involving the emission 
and re-absorption of the same photon by the same electron, and it is clear that 
such a transition leads to an infinite self-energy term in the solution. No such 
transition arises in the simple treatment of the propagation problem. Having 
excluded these roundabout transitions, it is obvious that, if two distinct transitions 
lead from an observable initial state to an observable final state, there is no reason 
to neglect one of them. 

The idea of “intermediate” states in the damping theory requires more 
detailed examination. The concept first arose in the treatment of the intermediate 
states occurring in Compton scattering. Those states have an energy which 
differs by a finite amount from the energy of the initial state for the scattering 
problem. Moreover, the energy spread of the final states tends to zero, as the 
enclosing “periodic” box becomes infinite; so the intermediate states cannot 
have the same energy as any of the final states. For scattering problems, therefore, 
the intermediate states are not observable; they have a mathematical meaning 
only. As a consequence, in the damping theory solution of scattering, the 


9 * This paper forms parts of a thesis presented to the University of Manchester for the Ph.D. 
degree. 


Damping Theory and the Propagation of Radiation | — 13 


possibility of the energy of the coupled system of particle and radiation being 
equal to the energy of an intermediate state is neglected. Further, in normalizing 
the total probability to unity, the possibility of the system being found in any 
intermediate state is neglected. 

The intermediate states (states in which one photon is present and both atoms. 
are unexcited) which may occur in the propagation problem are of quite different 
- nature from those which have just been discussed. Any of the intermediate 
states in the propagation problem may be observed; so they are no more intermed--: 
iate than they are final states. That they happen to be coupled to final states. 
of another kind is incidental. Further, it is clear that they have to be considered 
in normalizing the total probability; and it will appear from the detailed solution 
that they give the major contribution to the normalization equation. There is, 
therefore, no basis for treating these states by the same method as was used for 
the intermediate states in the scattering problem. 

It is shown that these considerations lead to a reasonable result for the pro- 
pagation problem. The method of discrete energy levels developed in an 
earlier paper (Hamilton 1947), to be referred to as I, is used, and the treatment 
of the simple atomic emitter is followed closely. The presence of the second 
atom causes a splitting of the energy levels of the interacting system, a splitting 
which is inversely proportional to the separation of the atoms. This splitting 
makes it possible for the second atom to remain in its ground state for the correct 
time after the first atom starts radiating. 

On the basis of Weisskopf and Wigner’s solution for the atomic emitter, 
Kikuchi (1930) has shown that the electromagnetic disturbance which is emitted 
travels with the velocity of light. Wentzel (1933) shows in a similar manner 
how the propagation problem can be treated. The present solution differs from 
these in that the reaction of the second atom on the first is considered, and that 
the solution obtained for the excitation probability of the second atom is a contin- 
uous function of the time. Both the above authors make substitutions for the 
excitation of the second atom which are discontinuous in time. 


§2. THE FUNDAMENTAL EQUATIONS 


The system composed of two identical atoms A and B and any radiation is 
enclosed in a large cubic ‘“‘periodic”’ box. It will be assumed that each atom 
has only one excited state, and states in which more than one photon occurs will 
be neglected. As was shown in I, the unperturbed states of such a system 
have a discrete energy spectrum; and it can be assumed that the distinct energy 
levels (which are highly degenerate) are equally spaced. 

The unperturbed states to be considered are: (0) atom A excited and no 
photon present, (1) atom B excited and no photon present, (i) both atoms in their 
normal state and a photon “i” present. 0, 1, i, will be used to label these states. 
The effect of the coupling between the atoms and the radiation field can be 
expressed by the equations (cf. I, eqn. (11)) 


(Ey -A)Ay + 2H Ai =9, 
(E;—A)Aj+ My Agt+HyA,=0, 7. eee (1), 
(E, —A)A, + 2H;A;=0. 
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‘The E denote the energies of the uncoupled states, the H are the usual matrix 
elements for emission and absorption of radiation, and the A are the amplitude 
factors describing the coupled states; A is the energy value of the coupled system. 
Different solutions of the set (1) will be denoted by the index p. 

Following I, any state of the coupled system can be described by the wave 
function “’, which can be written in the form ‘Y= &Xc,,‘V’" where the c, are constants 


depending on the initial conditions and We =A,“ exp {—i(A, —E,)t},. The 
index r runs over 0, 1 and all the i, while the %, are the wave functions of the 
uncoupled system. ‘’ can also be written 


es Gh a Sk eres (2) 


where a,(t) =A," exp {—7(A,,—E,)t} is time dependent. If the initial conditions 
-are oe 


a(O=1,  4,(0)=0, <a(0= 0p e= = eee (2’) 


then c,= A)". 

Due to the positions of atoms A and B being different (both atoms are regarded 
as being massive and immovable), the matrix elements H);, H,; are not identical. 
"They are related thus 

Hi Sexphi( peek )hijne =e ae (3) 


where r’ is the displacement vector of B relative to A and p, is the momentum 
of photon “i” (4 =c=1 throughout). 
Elimination of the A, from (1) gives 


Hal? | 4,» Hub 
{Hy-A+E ET bay + EEA, =0, 

os |Aail? Aulligg | 
{By ee eine eee yay teat | 


"Taking Ey =£;, these equations give 


| Foil?) Haj? — AoA} _ 


; | Hy)? { 
(EA) +(E, NEG EE TEVA) een, (4) 


where the bar denotes the complex conjugate. 

Using (3), any term of numerator of the last summation in (4) becomes 
| Ho;|?| Ho;|?. [1 — exp {—2(p;. r’)} exp {a(p,.r’)}]. . The summation = in (4) can 

i 
be split up into a summation over the discrete energy levels, denoted by ¥, and 
a summation over all the states belonging to a fixed value of E, denoted by 3 
. . . . . [#i] 

Consider, for simplicity, that the excited state of either atom is a p-state with 

the magnetic quantum number zero. Then 


Hy =const. | g(r)(z/r) exp {i(p;-x)}f(7) dV, 
with an obvious notation. Thus for dipole radiation 


Hy =ip;,..const. | g(r)f(r)(2"/r) dV, 
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where p,, is the z-component of p;._ Hence 
| Hul? =(pi.)°C, 
where C is independent of p,. 


Denoting the angle between the axis OZ and r’ by 6’, the angular summation 
of the term containing the exponential factor gives 
sin 


B exp { — i(p,.2")} Hol? =H. 4 + (3.cos?6’—1) Cr Hey Ones eh, 
[Fi] of 


a a3 


where | Hy;|?= & | Hy,|? is the summation of | Hy,|? over all angles, and « =| p;|r’. 
Ej 
If the zero of energy is taken as the energy of either of the atoms in its normal 
state, then | p;| =£;. 
It is convenient to write the angular summation in the form 


2 exp {— 1(p;- r’)} | Hoi]? =| Aol? {G,(Ejr’) sin(E;r’) + G,(E,r’) cos (Eyr')}, 
[4i] 


where 
|G (Er) =3 cos? 6’ /(E;r’) —3(3 cos? 0’ —1)/(E,r’)?; G.(E,r') =3(3 cos? & —1)/(E,r’)?. 
Then 
X & | Aoj|*| Hol? [1 — exp { —2(p;. r’)} exp {2(p;. 2’) }] 
[4] [43] Beppe tee 
=| Ao,|?. | Hoi? (1 —[G4(Z;7’) sin Ey’) + G,(E,7’) cos (E7’)] 


x [Giger ysin (G7) + G.(B 7’) cos(Eyr’)]}.. eee. (5) 


§3. THE ENERGY SUMMATION AND THE 
EIGENVALUE EQUATION 
Assuming that the matrix elements H; vary slowly with the energy £, in the 
important region (i.e. over the ordinary line width) it is possible to separate off 
energy summations in equation (4) of the type X1/(A—£;). Assuming that the 
E: 


distinct energy levels EF, are equally spaced in the vicinity of A such a summation 
gives 
x 1/(A—£;)=(a7/AE) cot (A—E,)n/AF}+o, (6) 
ij 


where AE is the spacing of the energy values, and £, is one of the energy values 
E, which lies close to A. o is a correction arising from the E, values which lie 

‘farfrom A. As the range of F; values centred about A over which the summation 
is made is increased, a is zero at first, and only becomes appreciable for ranges 
which are very many times greater than the l’ne width of one of the atoms. Over 
such large ranges the matrix elements will vary appreciably with £;, and (6) 
would have to be modified. In the extreme case of an infinite range of FE, values 
the term oa leads to the usual transverse self-energy of the bound electron. It 
is sufficient here to note that in the present problem o can be neglected for ranges 
of E, which are many multiples of the line width. Such a limitation will be 
used in what follows. Under similar conditions 


SI = a / =e / 
y sin (A E;)r sei 5 cos (A wr 


“SS a =(/AB) cot {(A -E,)n/AB}. 
Ej en EB; Wire SF 
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The variation of G,(Ejr’), G.(E,r’) with E, in the energy summations occurring 
in (4) can be neglected, provided the “periodic box”? is sufficiently large. These 
functions are therefore replaced by G,(Ar’), G,(Ar’) respectively. 
Using these summations, (4) with the aid of (5) reduces to a quadratic equation 
for (7/AE) . | Ho,|2. cot {(A — E,)7/AE} whose roots are 
A-—E, ¥ {G,(Ar’)sin (Ar’) — G,(Ar’) cos (Az’)} . (1/2) 
1 + {G,(Ar’) sin (Ar’) + G,(Ar’) cos (Ar’)} ; 


(I'/2) cot {(A —E,)r/AE} = 


where the upper or lower signs are to be taken together, and ['/2 =(7/AE). | Hy\|*. 
I is the ordinary line breadth of either atom (cf. I, §5). According to (7) there 
will be two roots A lying in the energy range AE between each two adjacent 
distinct values of E;. As r’-> 00 equation (7) becomes 
(1/2) cot {(A —E,)7/AE} =A —E,, BARS 

which is the eigenvalue equation for the photon emission of one atom. By the 
definition of the functions G it follows that 

G(Ar')~(1/Ar’) O11), Ga(Ar’)~(1/Ar’? O(1), 
provided Ar’>1. O(1) denotes a constant of the order of magnitude unity. 
Thus the roots of (7) differ from the roots of (7’) by amounts of the order of 
AE/(Ar'), provided Ar’>1. As 1/Ar’=A/r’, where A is the wavelength of the 
resonant radiation, this condition can be written as 7’ >A. 


§4. NORMALIZATION 
The condition that the total probability of the system being in all configurations 
is unity, gives 


[Age Ae Age = tA) eee oe (8) 
for all values of the root index ». Substituting from (1) gives 
HH; We erle 1 
A,’ =4>—*—% x ae 
tn Rg aot Og AN) Oe, 


Making these summations by the method of §3, and using (7) to evaluate 
cot {(A — E,)z/AE} where necessary, it appears that 


Ait= ots Aa, ene! wisi (9) 
where the + sign is to be taken according as the upper or lower set of signs is 
used in (7). Further, from (1), (3) and (9) 

%| Ap? =| Ay? | Hyl?.1 exp (—i(p,. v)}P/(A—B) 


=2| Ag |?.| Hy\|?. (7/AE)? . cosec? {(A — E.)r/AE} 
x {1 + [G,(A7’) sin (Ar’) + G,(Ar’)cos(Ar')]}, ... (10) 
provided 7’< I, the linear dimension of the enclosing box. 

Now it is clear, by using the relation ['/2 =(z/AE) . | Ho;|?, that in (8) the terms 
| Ay"|?, | A,"|? are unimportant, as 7/AE is very large. The normalization is 
therefore determined by the photon states. From (10) and (7) it follows that 
1 | Hol? {1 + [G, sin (Ar’) + G, cos (Ar’)]} 
2° {A—E, ¥ (IT'/2)[G, sin (Ar’) — G, cos (Ar’)|}? 

+ (T/2)2{1 + [G, sin (Ar’) + G, cos (A7’)] 2 
heer (11) 


| Ag" |? = 
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So DHE EXCITATION OF LBRE ATOMS 


It follows from (2) that the probability of the occupation of the state ys is: 
given by | a|?, where ap is related to the Ay“ through the equation 


age fae Ao" exp —1u(A,,—Ep)t}. 


For the initial conditions (25) o Hy so 
a(t) exp { —iE ot} = X[ Ap"? exp { —iA“2}, 
7 


where the summation runs over all the roots of (1). Using (11) it is seen that 
the summation over the pair of A, values which lie in each energy interval AE 
just cancels out the terms of order A/7’, leaving a large term, plus terms of order 
(A/r’)?. (The slight difference between the two values of A, give corrections 
which vanish as Loo). Omitting the terms of order (A/r’)? the sum of 
such pairs of | A,"|? values is (|A "|?)’ =| Hy,|2/{(A —2Z,)?+02/4};. (| Ag“ |?) is 
identical with the term | A,“|? which arises in the single atom problem (cf. I). 
Replacing the summation in (12) by an integration it is readily seen that 
| ao(t)|? =exp (—T?) as for the single atom. 

The behaviour of the probability function | a,|? for the excitation of atom B 
is quite different. ‘This function is given by 


a,(t) exp (—7E,t) = BAA Seep Ct Agh)s GL” Sea (13) 


Relation (9) causes a cancellation of a largest term in | Ao“|?, on adding the 
contribution from the pair of values within a range AF. ‘The remaining terms 
are smaller than (| Ay“|?)’ by a factor of the order A/7’. 

Equation (13) can be written in the form 


a, exp(—iE,t) = =| A,*+|?exp(—zA,,t)— 2 | Ap |2exp(—7A,_2), ...... (14) 
BM u— 


where the +, — indicate the two types of roots corresponding to A,“= + Ay. 
Also, (11) can be written 

| Agi}? = £ Feo ey earl PE (YPH CGHESLENIST ca 
. (A,,— £,—10'/2) + (1'/2)(G, —1G,) exp (—2A,7’) 


1 
~ (A, — FE, +11/2) + (1'/2)(G, +7G,) exp aral 
It can be assumed that A,,, =A, =A, as the error introduced thereby is negligible. 
Replacing the summation in (14) by an integration gives 


ass =a eee Gece aC) 
pot) 5 \5 (N—E, +102)? (1 /2)°(G, +7) exp (Ar) 


(G, —iG,) exp {—zA(r’ + t)} | 
(A — £, —11'/2)? — (1'/2)?(G, —17@,)? exp (— 2k) 


is 


Equation (15) is exact, provided the dimensions of the “periodic box” are 
great enough; and this equation could be used to calculate a, to any power of 
(A/r’). Here it will be considered sufficient to find a, to the lowest power of 
(A/r’). Thus Gzisneglected. Further, the variation of G,(A7’) with A is neglected 
for the range of values of A which is important; and from (15) it is clear that this 
range is of the order of I’. 
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With these assumptions the integrals in (15) can be evaluated by contour 
integration. The circuit consists of the real axis, together with half of the circle 
at infinity. On the circle at infinity the second integrand takes the form 


exp {y(r' +1)}/(R? +Cexp (2yr’)), where A =x+1y, (aN ve. and |C| is a constant. | 


On the half circle y <0 this integrand will give no contribution, provided t>—r’. 
Moreover, the poles of the second integrand all lie on the half plane y>0. To 
see this, substitute A’ =A —E,—7I'/2. The poles are given by 

A’ = +(0'/2)G, exp {((L'/2)r' —i(E, +A’)r}. 
If A’ =x’ —1y’, this gives 

V(x’? +y"*) =(T'/2)| Gy] exp (1/2 — yy}. 
As |G,|<1, y’>I/2 is impossible, and hence the result. Thus the second 
integrand gives no contribution to a,, provided t>—r’. Actually it represents 
the “incoming” wave and plays no part in the solution for ¢>0. 

It is readily seen that the roots of the denominator of the first integrand all 
lie on the half plane y<0. On the infinite circle the form of this integrand is 
exp {—y(r’ —t)}/{R?+Cexp(—2yr’)}. When t<7’ this gives no contribution 
when integrated over the upper half circle, so a,(¢) =0 if t<7’. 

For t>r’ the value of the integral can be found by noting that, as there are 
no roots of the denominator on the upper half plane, the integral is equal to 

a Gy FIN ot AOSTA (a8) eee a ieee 
ce (A—£,4+70/2)? —(T/2)?G,? exp (2zAr’) 
where C is the upper half circle at infinity. On C the second term in the 
denominator can be neglected, giving — Gi] dA [exp {2A(r’ —2)}]/(A —&, +71 /2). 


The path C can be closed by including the lower half circle at infinity which gives 


no contribution when t>7’. ‘Thus the integral becomes 
— Gy | aA [exp ftA(r’ —1)}]/(A —E, +i /2), 
where O is any contour surrounding the point (£, —7I'/2). Hence 
a,(t) exp (—2E,2) =1G;.(I'/2)(t—7’) exp { -—T\(¢—1’)/2} exp (iE,r’) — (¢>17’), 


and as | G,|?=K(A/r’)?, where K is of the order of magnitude unity, the final 
result is 


| a,(z) |? =0 t<r’ 
|a,(0) 2 = K jr {PC 1')2}2exp{-L(e—r)} | 
Thus the atom B remains in its normal state till time 7’, after which it is excited 
for a time of the order 1/[. The maximum probability of excitation is of the 
form Kar’), which is in agreement with the usual ideas on electromagnetic 
propagation. 

The resonance phenomenon will arise from the re-emission of radiation by 
the atom B, and its absorption by A. The excitation amplitude of A will be of 
the order (A/r’)?, and therefore does not occur in the approximation used in this 
section. 
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ABSTRACT. The cross-section for Bremsstrahlung in neutron-proton collisions is 
calculated with a central exchange-force model and for energies up to 20 Mev. Its value 
is found to lie between 10-8 and 10-2 cm?. This probably is too small to be measured 
with experimental techniques at present available. 


SiN ERO DUC LTON 


HE aim of this paper is to estimate the magnitude of the cross-section for 
| Bremsstrahlung in n—p (neutron—proton) collisions and its sensitivity to 
variation of the hypotheses concerning the n—p interaction. Neutron 
and proton will be treated here as point particles interacting by means of a central 
force of the exchange type. ‘Their interaction with radiation will be supposed to 
conform to the usual theory for particles in a radiation field. We thus omit 
refinements which have been included by Pais (1943) in his treatment of the photo- 
disintegration of the deuteron. 
This simple model is inadequate when the incident neutron energies are high 
(of the order of the meson rest-energy). Our discussion will thus be restricted to a 
lower energy-range. The relevant energies will then, in general, not be large 
compared with a suitable average of the potential. In these circumstances, use of 
the Born approximation may result in appreciable error. For this reason, exact 
continuum wave functions of the n—p system will be used to calculate the required 
cross-sections. 
§2. WAVE FUNCTIONS 


Referred to the centre of gravity frame, (C-system), the wave equation for the 

internal motion of the n—p system is 

[V2 +(M/h?)(E —V)]b(r; my, m,) =0. pee (zel:) 
Here r=(x, y, 2) is the position-vector of the neutron relative to the proton, m, 
and mz, are spin-variables of the particles, MM the mass of a nucleon and 2£ is the 
kinetic energy of the incident neutron in the laboratory system. ‘The interaction 
operator V will be chosen to represent a mixture of central exchange forces of 
Majorana and Heisenberg type (cf. Bethe and Bacher 1936, §§11, 12, 13). 

The wave functions will be required to represent “incoming” or “‘outgoing”’ 
particles which have, at large distances, a definite momentum #k. Correspond- 
ingly we determine the following types of solution of (2.1): 

(a) To represent incoming particles: solutions which have, for r— o, the 

asymptotic form 
[exp {7(k.r)}+ diverging spherical wave] x (spin function) 
(b) To represent outgoing particles we choose, for reasons of symmetry, those 
solutions which have the asymptotic form 
[exp {7(k.r)}+ converging spherical wave] x (spin function). 
2-2 
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The eigenstates can be classified into triplet states with wave functions 
w= UME) ty a0, el oy eee (2.2) 
and singlet states with wave functions : 
eas Ohi = ep Cea ee Bg) wen ir (233) 


3, (m=0, +1) and 1y are the usual normalized symmetric and antisymmetric 
spin functions of the two-particle problem. The wave equations for triplet and 
for singlet states are independent. Also, the space part U(r) of the triplet wave 
function (2.2) is independent of m. On separating out the spin-dependence, the 
functions U® and U® satisfy equations which may be written together in the form 


(AIM +E)\U@)=V nO 2)s ye eee eee (2.4) 
If M(r) and H(r) describe the dependence on distance of the Majorana and 
Heisenberg potentials, respectively, then 
M(r)+H(r)_ (triplet states) | 
V(r) = : : 
M(r)—H(r) (singlet states) 


Each wave function will be characterized by an associated wave vector k of 
length k=(ME/h?). This dependence on the parameter k will be indicated by 
the notation U(r; k). Introducing spherical polar coordinates (r, 0, 6) with polar 
axis in the direction k, those wave functions which satisfy conditions (a) or (6) 
above can be cua in the form: 


U.(r;-k)= males 1)’exp { +2,(k)} —-— 


ts 


ude sh) P(cos@) —...... (2.6) 


Pet conthg! dpe ato tame enter ely ‘lee al 
rs = nee eel 2 


(cf. Mott and Massey 1933), The P(t) are hanes: Legendre polynomials 
and u(r; k) is that solution of 


| 
du + [a a -(-1 7] a eEt bape (2.7) 


dr® 
which vanishes for r=0 and which has, for r>0o0, the asymptotic form 4 
u(r; k)~sin[kr—fatn(k)) cae (2.8) i 


U., and U_ correspond respectively to the cases with diverging and conve: sue 
Sue waves. 

The incident beam will be assumed to be unpolarized and to consist of particles 
with wave vector ky (in the C-system). There are then four independent initial 
states of equal weight : 

UEXE; ky) eX (m=0, + 1) 
Pit; Ky)= | e SS (2.9) 
UN"(t; Ko) *x J 


U® and U) are given by (2.6) with k=ky. Similarly, there are four independent 
wave functions corresponding to final states with wave vector k: 


; Bog OSes ie ‘Xm (m=0 Pp 1) 
bir (r; k ra UM; é 


U®) and U™ are given by (2.6) with z as a axis. (Since V contains no spin- 
orbit coupling terms, the spin functions in both (2.9) and (2.10) may be referred. 
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to the same coordinate system with ky as z-axis.) We shall denote coordinates 
in asystem with ky as z-axis by (x,y,z) or (7,6,4) and polar coordinates in a 
system with k as polar axis by (7, 6’, ¢’). 

Each of the functions in (2.9), (2.10) is normalized per unit volume, i.e. 


col ele eee 


Ui CO: v 


§3. MATRIX ELEMENTS 


We shall consider processes in which the n—p system (referred to the C-frame), 
makes a radiative transition from one of the states characterized by wave vector ky 
(energy E)) to states characterized by wave vector k (energy FE). The circular 
frequency w of the photon emitted is then given by 
ho=Hy—B=(hiM\Reak). a. (321) 
Strictly speaking, the C-systems appropriate to initial and to final states are not 
identical. ‘The difference—due to the momentum carried off by the photon—is, 
however, very small and will be neglected. 

Two types of process will be taken into account. These are (Bethe and 
Bacher 1936, §§16, 17): (a) transitions induced by the electric dipole der; and 
(2) transitions induced by the magnetic dipole wpT=po(u, 6) + HrGn), Sp Sy are 
the spin operators and yu, 4, the magnitudes of the magnetic moments in units of 
a nuclear, magneton py =eh/2Mc. 

Correspondingly, we shall have to evaluate matrix elements of the following 
two main types: 


MC =heM=he | Ph (rs k)rpn(t; ky)dr, wees (3.2) 
N= poN =o | oat(rs KT Y(t; Kolar. eee (3.3) 


The required cross-sections are to be obtained by calculating the eight partial 
cross-sections with each of the states (2.9) as initial state and then adding them 
with weight-factor + for each. | 

Considering first transitions of type (a), the only non-vanjshing matrix elements 
are: 

(i) three equal elements : 


M®) — | UP(£; ky)rUO*(r; k)dr, Be ens.) 


corresponding to transitions U®.®y,,>U®). *y,,, (m=0, +1), and 
(ii) one element : 
M® = | UM(r: kyrU Or; k)dt, cae (3.5) 

corresponding to the transition U(P!y >U® Ly, 

In the sum for the composite cross-section, (3.4) contributes with weight ? 
and (3.5) with weight 4. 

Let (A, ») be the angular coordinates of k in the coordinate system with kp as 
polar axis. For conciseness we write 


P=1P®, APH), F=3 orl) sees (3.6) 


and 


Gi n= ale RNIN TPE C0 at alae ales (On) 
ott Rok 0 
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In (3.7) G®, and G®, are to be distinguished according as triplet or singlet wave 
functions and phases are involved. 

As it stands, the integral in (3.7) is divergent. However, it is summable, e.g. 
by inserting a convergence factor e~ and then proceeding to the limit b— 0 after 
the integration has been performed. _ It is in this sense that (3.7) and the analogous 
divergent integrals in (3.13) have to be interpreted here. 

When the wave functions in the form (2.6) are substituted in (3.4) and (3.5) 
the integrals over angles are easily evaluated by using the spherical harmonic 
addition theorem. We find (writing t=cos A) 


M,=4r1 > (1+ 1){Gyy, exp Fras +} PD) — Gyisa exp H+ Ma) PO}, 


1=0 


M,, +1M, = —47 exp (+t) cr (Giya,1€8P WHi44 +m) P(t) 


+ Gyr exp (q+ Mad PiiO}. «+. (a0) 


Turning now to the photomagnetic transitions, the et surviving elements 


Nk); k) are of the forms: 


NC®-» corresponding to transitions triplet— singlet, 


NC*® corresponding to transitions singlet — triplet. 


Writing 
Abd = (UE; U(r hdr a © [Gp lee(set oor iss in eee (3.10) 
we have 
> {INS D2 a DA—2 YF | ND 2—%y —p)?| AGIA ...0.0, 
ere NL oe detailed (Hp — Bn)” | A%?)| (3.11) : 


The cross-sections obtained from (3.11) with (7,7) =(3, 1) and (1, 3) each contribute 
with weight + to the final cross-section. 


Sbetinine wave functions of the form (2.6) in (3.10) we find analogously to 
(3.3) (9.9)2 


AG) =4q x ACESS: Dexp {i {un 7H +7) P(t), are iagat ee (3.12) 


with 
ij ripe i j 
KG = aed UNS Ry (wk) area eee (3:13) 


§4. CROSS-SECTION FORMULAE 


Each wave function (2.9) for an initial state corresponds to an incident current 
of 2hko|/M neutrons/cm?/second. Each final wave function (2.10) is normalized 
per unit volume; there are then (Mk/167°h?)dE dQ. final states of each type 
corresponding to motion within dQ=sinAdAdw and with energy in the range 
(E,E+dE). The factor 


Me? ; 
aan -) dBA cari) Une 5 Seer (4.1) 


is thus required to translate transition probabilities calculated with the wave 
functions (2.9), (2.10) into differential cross-sections. 
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Let o(£y; £,A,«)sinAsinadAd« be the cross-section for a process in which 
(i) the n—-p system makes a radiative transition from a state (ky, Ey) to a state with 
energy in (LE, +dE) and direction of motion making an angle between X and 
A+dd with ky; (ii) a photon with circular frequency w =(Ey —E)/h is emitted in a 
direction making an angle between « and «+da with ky. Then 


M? (e} 


3 3 
oly: IOI ais 32h z=) pall Me +1M,)(1 + cos?) +2|.@,|? sin? «}, 


rere AC is to be interpreted as AC or AC in the electric dipole case.and as 
NN» or A(* in the magnetic dipole case. The corresponding cross-sections 
will pe distinguished by affixing the same indices to o. 

Further cross-sections of interest are derived from (4.2) by integration. 
They may be written in an obvious notation as o(Ey; E,A), o(Ey; E; «), o(Ey; E) 
and o(£,). In particular 


WP (E\* : 
alls; EZ iam E) ), 2. |.&@,.|2sinAda. Sens, aNore (4.3) 
From (4.3) with (3.8), (3.9) and with (3.11), (3.12) we have 
; 2e?M7w? ( E\? 
oes B)= Sree (ze) ECHR HCI (4.4) 
Os § Ly — Uy)? M2 3 E ye A me 
ot )(Ey; E) = pole (=) Syl PBSC nes (4.5) 
o/ 1=0 


Soe SORT WANG EH TEORCES 


For further progress it is necessary to introduce explicit assumptions for the 
interaction function V(r) in (2.4), (2.7). In this section we assume only that the 
forces are of short range, so that we may put effectively 


VACA) Rs 2 7 ED OR a (5:1) 


In the range r>a, the required solution of (2.7) can be expressed in terms of 
Hankel functions: 


u(r; k)=4 (7 “y {exp (tm) H{),(kr) +exp(—in) HPP (Rr)}. wee eee (5.2) 


The phase 7,(k) is determined by fitting the interior and exterior solutions smoothly 
at r=a; it depends on the detailed form of V(r) in (0,a@). (5.2) can also be repre- 
sented in the form 


u(r; k)=F(kr) sin (kr +7) +G (kr) cos(kr+).  ..-0- (533) 


F(x) and G,(x) are polynomials in x1, defined through the Bessel functions by 
($27)! J,,,(x) = F(x) sin x + Gi(x) cos x. 

Now, the hypothesis of short-range forces implies that, given k, there exists 
an integer /,(k) such that for all /> Jp, (i) m,(k) does not differ appreciably from zero, 
and (ii) the u(r; k) do not differ appreciably from the unmodified plane-wave 
components (zkr/2)'J,,,(kr) in the whole range (0, 0). It follows that, for 
1,n >1)(Ro), the G,,, in (3.7) can effectively be taken equal to 0. Similar consider- 
ations apply to the Kj. 


24 M. Krook 


We now introduce the approximation of replacing the integrations over (0, «) 
in (3.7) and (3.13) by integrations over (a,). This approximation is valid only 
for energies which are not too high; this has been checked for the case of a 
rectangular potential (cf. §7). For energies much larger than those considered 
here the contribution from (0, a) is no longer negligible; the calculated cross- 
sections would then depend more intimately on the detailed form of V(r) than in 
the low-energy case. 

The integrals for G, ;,, and K, are now easily Malate by elementary methods. 
(The evaluation for higher / can be simplified by using for the u, certain recurrence 
formulae analogous to, and easily derivable from, those for the Hankel functions.) 
In particular we have 


5 E v\t $ 
(ko — k?)?Go 1 = (2 = 3) C5, +2 (=) SoC — Ro (1 a z) | (=) foC1 + 8 | , 
(5.4) 


eee eee 


: pe E ENS 8 E - Eee 
(kj —R7)G4, 0 = (= 3) C951 +2 (=) SoCal kya (I = rz) Ee + (3) sis , 


eee (5.5) 
where, for conciseness, we have written (cf. (3.6)) 

=s=sin(ka+7q®), =? =sin(ka+n®), etc. ...... (5.6) 

Also 

Ql<¢ 3 3 S « 
— RPK O = (1 ze zm) | (3) asp) He | PES) (5.7) 
2 ENE AE \ Ae 

(2 RyPeK 9 = (1 2 z) | (3) ay — sep |. wee (5.8) 


q 
§6. FORCES OF “ZERO RANGE” i 
We consider the simple limiting case where the range a>0. Here ») 0 for f 
all/>0. Let «; and —e, be the binding energies of the deuteron for the ground i 
state and for the virtual singlet level respectively. Then : 
ka +1) $n + arc tan (e,/E)!; ka +P $a —arctan(e,/E)*. ...... (6.1) i 
In this case G, 141, Gj41,, and the K, are zero for all J40 and 
E \t Eo | 
he? — p2)2 Pa , 2 . 7 
(RC, 2(—) , (BURG, = 2(- =)" Snel (6.2) 
e being taken as ¢, for the G® and as «, for the GY. - Also 
~— 2 
Rk? — R2)3/2K (49) = (et 4 3 k Ey-# 
( 0 ) 0 (€3} +€}) (e; + Eye; + £) es . Melee ale (6.3) 
The various cross-sections are then given by 
2:13 x10- /E\t/( 4F 4E 
o(Ey; E)= red pan ails o | 
Sa 40 E, 25 (= \e torhh? So e.cic (6.4) 
o@\(Ey; £)=1-63 x 10-3 G i fi te eae 
Lo 5 (<} +e4) (ERE) oe (6.5) 
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Table 1 lists total cross-sections (integrated over £) for the various types of 
process, and total averaged cross-sections o(£), for a few values of Ey. The integral 


for o,(Z) diverges logarithmically at E,; we have taken 
SPE, 
o{fy) = | » Ey; E) dE. 


‘The calculation has been made with «,=2-17, «;=0-05 Mev. in (6.4), (6.5) and 
with f2=0-95. 


‘Table 1. Total Cross-sections in Units of 10-2°cm?2 for Radiative Scattering by 
Protons of Neutrons with Energy 2E, Mev. 


Ey (Mev.) a op mag o(E) 
1 119 3°26 0:02 ey, 
4 2°64 4-30 0-05 aig | 
6 3:07 4°50 0-06 3:5 
10 3°56 4-67 0-08 oa 


§7. FINITE RANGE 

The cross-sections were also calculated for energy Ey =10 Mev. using a rect- 
angular potential model with range a=2-5 x 10-%cm., (V),—, =25 Mev. (triplet 
states), 15 Mev. (singlet states). By far the largest contribution is again due to the 
electric dipole transitions /=0—/=1. 'The approximation of replacing inte- 
grations over (0, 00) by ones over (a, 0) (cf. §5) was checked for this case. Inno 
case did the neglected contribution from (0, a) amount to more than 10% of the 
total; in general it is considerably less, as also is its averaged contribution. In 
Table 2 are listed the cross-sections for a number of values of E (column 2). For 
comparison, o(10, £) in the limiting case a=0 are given in column 3. 


‘Table 2. Partial Cross-sections in Units of 10-?° cm? for Radiative Scattering by 
Protons of Neutrons with Initial Energy 20 mev. 


_ E(mev.) o(10, E) o(10, E), (a=0) o(10, E), (max) 
1 1-7 1:3 26 
5 2-7 2:3 18 
al 47 4-2 29 
9 15-0. 13-8 


We have noted that the major contributions to the cross-section, for initial 
energies Hy) up to about 10 Mev., arises from the terms in Gy , and G,,) ofo,. For 
such energies Gy , and G, , are closely approximated to by (5.4) and (5.5) irrespec- 
tive of the detailed form of V(r) in (0, a). We thus obtain an upper limit to these 
coefficients, and so to the calculated cross-sections, by replacing all sine and cosine 


factors by unity. ‘Then 


E EY} E\[/E,\! 
(3-H) Gai <5 -3|+2(F] +ha(1~=)| (3) +1], 


p3/+2() tee(1-z) [(Z) #1] 
B _3)2(2)' +4ee(i—Z)[(2) +1] 


Column 4 of Table 2 gives maximum possible values of the cross-section obtained 
by using (7.1) with EZ) =10 Mev. 


(kj — k?)?| Gio! < 
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The cross-sections for Bremsstrahlung on the central exchange force model 
thus come out very small in the energy range considered; in n—p collisions there 
would be about 10,000 elastic scattering processes for each such inelastic process. 
The relative probability for a given direction of the proton has also been investi- 
gated and found to be of the same order of magnitude. It thus appears that it 
would not be possible to detect and measure such radiative collisions with 
experimental techniques at present available. (In view of the smallness of the 
cross-sections it would be pedantic to discuss here the question of angular distri- 
bution of the particles and photons.) For much higher energies the cross-section 
for Bremsstrahlung would increase appreciably (at least on this model), and might 
be capable of measurement. The model used in this paper would, however, be 
inadequate for treating the problem at high energies with any exactness. 
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ABSTRACT. ‘The resonance transfer of energy in crystals of organic compounds when. 
excited by visible or ultra-violet light has been examined by measuring the fluorescence 
of the systems naphthalene—anthracene, anthracene—naphthacene, naphthacene—pentacene, 
anthracene—acridine, anthracene—phenazine, and’ anthracene—di- and hexabromobenzenes. 
The conditions for ‘‘ exciton’ trapping and the relations of resonance transfer to fluores- 
cence quenching are discussed. 


SG IBN AEICONONOIE MIKO! 


HE transfer of excitation energy by resonance from molecule to molecule in 
| organic crystals is most strikingly shown by the anthracene-naphthacene 
system (Bowen 1945, Bowen and Mikiewicz 1947). Very small quantities of 
naphthacene in solid solution in anthracene change the fluorescence from blue to 
green, the excitation energy travelling by exchange until it is trapped at a 
naphthacene molecule. ‘This effect disappears when the crystals are dissolved 
in benzene. ‘The resonance exchange is doubtless assisted by the fact that the flat 
molecules of anthracene lie nearly parallel in the crystal (Robertson 1932); in the 
solid solutions the naphthacene molecules are again parallel to them. 


: 
| 
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For the unambiguous demonstration of the phenomenon it is necessary to use 
systems where the extinction coefficients show (a) that the exciting light is absorbed 
by the major constituent, and (d) that the fluorescence of the major constituent is not 
reduced through secondary absorption by the added substance. For ‘“‘exciton’” 
trapping to occur the added substance should have a possible transition between 
energy states associated with a lower energy quantum than a similar transition in 
the major substance, i.e., should absorb at somewhat longer wavelengths. This. 
implies that the fluorescence band of the latter will almost coincide with the 
absorption band of the former, but condition (b) will hold if the added substance 
acts at sufficiently small concentrations. Extinction coefficient measurements are 
usually wanting for the crystalline state. The crystals are optically highly aniso- 
tropic. From the results of Kortum and Finckh (1942) on anthracene, however, 
it appears that there is little difference between the averaged values for crystals and 
solutions, except for a shift of about 200. of the band to the red in the case of 
crystals. - Solution values have therefore been used after slight correction. 

A very high degree of purity is necessary in this work, and the materials used 
were subjected to careful purification and checked by absorption or fluorescence: 
examination. Mixtures were prepared by weighing the components into a small: 
tube, evacuating, and fusing together. The crystals were ground up and spread in 
a thin layer on the flat outer side of the inner part of a specially designed metal! 
Dewar vessel with a silica window for excitation and observation. The fluores- 
cence was excited by a concentrated beam from a small high-pressure mercury arc. 
with filters, and measured through other filters by a photo-electron multiplier unit.. 

The excitation energy in the crystal may undergo a variety of processes, which 
may be formally represented by the following scheme, in which the k’s are velocity 
constants, the square brackets indicate concentrations and A*, B* =molecules of 
A, B with excitation energy. 


Process Relative Rate Description of Process 
A +hv=A* 1 Excitation by light. 
Lakes ee a, k,[A*] Fluorescence of major constituent. 
PA — Al R,k[ A*] Self-quenching of major constituent. 
EES 8) seid) a Als Exciton transfer. 
A*+B =B*F+A k,k;[A*][B] Exciton transfer to added substance. 
JB} =B -+hr, k,{B*] ' Fluorescence of added substance. 
Bs ss k,k;| B¥] 
Be B= Bb 4B Ryko B*|[B] Self-quenching of added substance. 
AUS) 1B) Ss lB) =e zl k,k-[A*][B] Quenching by added substance. 


The quantum efficiencies of fluorescence are then for A: 


1/{1+k,+k,[B] +k,[B]} ; 
and for B: ' 
k[B]/{1 +k, +k,[B] + k[B]}{1 +s + Re BI}. 


SE ReSLI TS 
The following systems were examined and results are shown graphically in 
Figure 1. Ordinates are fluorescence efficiencies (quantum yields) and abscissae- 
the concentrations of the minor constituent expressed logarithmically (e.g. 
_7 indicates 10-7 gm/gm.) except in the case of the anthracene-bromobenzene 
systems, where they are given as gm. per cent. The values of molar extinction. 


28 E. #. Bowen, E. Mikiewicz and F. W. Smith 


coefficients « given below for comparison relate to solutions of the hydrocarbons 
with a slight correction for band shift in the solid substances. 


NAPHTHALENE —ANTHRACENE ANTHRACENE = peasaliieeces 
10 


Figure 1. Abscissae: concentration of minor constituent expressed logarith- 
mically (for anthracene—bromobenzene, expressed as a percentage). 
Ordinates: fluorescence efficiencies. 


ye A ING A SN 
Naphthalene | Se | with anthracene | | { |. 
7 Ve 
Exciting light: 2900-3200 ., filtered by solutions of nickel, cobalt, 


and copper sulphates. 
Absorption: naphthalene 2500-3200a.,\ loge =3-5 at 3000. 
anthracene 3200-4100a., loge =2-4 at 3000, and 3-5 at 3650. 
Fluorescence: naphthalene 3130-3720 a., estimated by amount of fluorescence 
cut off by an Ilford 805 filter. 
anthracene 4020-5020a., estimated from total fluorescence less 
that removed by an [ford 805 filter. 
Crystalline naphthalene has an efficiency of about 0-7 at ordinary temperature, 
-but in the anthracene solid solutions the anthracene fluorescence rises to an 
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efficiency of 0-9, equal to that of pure anthracene, at a concentration of 10%. 
At —180°c. the efficienceis are 0-95 and 1-0 respectively ; k,=0°3, k,=5. x 104 at 
20°c., kp =0:05, k,; =1-3 x 104 at —180°c. The other constants are zero, 


NY aS PAN GN a 
Anthracene | | | | with naphthacene | | ho ai 
SNS IAN NN 
Exciting light: 3650a., filtered by Chance’s glass OX 1. 


Absorption: anthracene 3200-4100a., loge =3-5 at 3650. 
naphthacene 3600-5000 4., loge =2-4 at 3650, and 4-0 at 4900. 
Fluorescence: anthracene 4020-5020a., estimated through Ilford filter 601. 
naphthacene 4900-5750 a., estimated through Ilford filter 404. 


The falling off of the naphthacene fluorescence at concentrations above 
about 0-1°% is attributed to self-quenching effects. The limit of anthracene—_ 
naphthacene solid solution formation is unknown, but appears to bé much greater 
than 0-1°%, so that it cannot be ascribed to the separation of a different solid phase. 
Pea Uleke= 2 <10° ko=2:5 x 10%at20°C.,.k-=0) kg = 0x 10%, ke = 107 at — 180° ; 
the other constants are zero. | 
oo eS os RENT BS 


le enh eacriaine = {7 ' | ae 
‘SEF SZ ONIN 
Exciting light: 3650 a. 
Absorption: anthracene 3200-4100a4., loge =3-5 at 3650. 
acridine 3000-4300 a., loge =3-6 at 3650. 


Fluorescence: anthracene 4020-5020a.; acridine feeble. 


Anthracene 


Since the solid solutions containing acridine show practically only anthracene 
fluorescence the acridine traps and quenches the exciton. k,=0-1, k,=5 x 10? at 
20°c.; the other constants are zero. 


Se FR 
Anthracene ( a | with phenazine | ( = | 


Soe SwWON OG 

_ Exciting light: 3650 a. 

Absorption: anthracene 3200-4100a., loge =3-5 at 3650. 
phenazine 3100-5000a., loge =4-0 at 3650. 

Fluorescence: anthracene 4020-5020a.; phenazine none. 


As acridine, phenazine acts as an exciton quencher. k,=0-1, k,=7 x10! at 
20° c., k, =0, k, =3 x 104 at —180°c.; the other constants are zero. 

The naphthalene—quinolene system is closely similar to this one except that 
its bands lie at shorter wavelengths. 


VOSS OES ZORA XS NH 
brapringcene. | 8 “wih pentacene* “||. |. fe) 

Pee SS NG A NS 
Exciting light: 4358a., filtered by solutions of cuprammonium 


sulphate and sodium nitrite. 
Absorption: naphthacene 3600-51004., loge =3-5 at 4358. 
pentacene 4300-6100a4., loge =2-6 at 4358 and 4-3 at 6000. 
‘Fluorescence: naphthacene 5000-5900 a.; pentacene beyond 6000. 


Owing to the feeble fluorescence of this system observations could be made only 
at liquid air temperatures, and even here is was possible only to make rough visual 
estimations. The absolute values of the fluorescence efficiencies given are very 
rough approximations: k, =2 x 10?, k; = 108, k; =3 x 10 approx. at —180°c. 
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Anthracene with p-dibromobenzene and hexabromobenzene. ; 

Exciting light 3650 a. 

Though it forms solid solutions with anthracene up to 30% concentration, 
p-dibromobenzene exerts little quenching effect. Hexabromobenzene behaves 
similarly. These substances have absorption bands at shorter wavelengths than 
3650. and appear unable to act as exciton traps, though they act as quenchers for 
anthracene fluorescence in liquid solutions (Bowen and Norton 1939). In the 
latter state van der Waals associations of anthracene molecules with bromo- 
compounds appear to be responsible for the quenching (Bowen, Barnes and 
Holliday 1947). ; 

In addition to the above systems the three-component system anthracene 
with small quantities of both naphthacene and phenazine was examined. ‘The 
phenazine diminishes the observed fluorescence coming from both the anthracene 
and the naphthacene. The two-component system constants given above were 
found to reproduce the fluorescences observed through Ilford filters 601 and 404 
(blue anthracene and green naphthacene emissions) within the limits of error. 
‘The agreement between observed and calculated results in all the above instances 
is only approximate, due to the inevitable inaccuracies of fluorescence measure- 
ments from crystals. Large errors occur through variations of light diffusion 
losses in different preparations. ‘The solid solutions also are far from uniform in 
composition. When a solid solution of naphthacene in anthracene is prepared, 
for example, either by fusion or by crystallization from a solvent, the first crystals 
deposited have a higher naphthacene concentration than the liquid, and as 
deposition proceeds a series of solid solutions of different compositions are formed 
successively. The only way out of this difficulty would be to prepare specimens 
by crystallizing a small fraction only from a liquid and to find its composition by 
subsequent analysis. This would require large quantities of highly purified 
materials which were beyond our resources. 


Sop ADISKS Wiss) OUIN| 


A quantity of great interest in relation to these systems is the “life”? of an 
exciton process. The mean lives of anthracene fluorescence in the solid and 
dissolved states are not accurately known. Wood (1921) found a value of about 
10-° sec. for crystalline anthracene, but recent unpublished work indicates that it is 
really very much smaller. The average life for dissolved anthracene, estimated by 
the polarization method in glycerol, appears to be 10’ sec., but this is probably too 
long because the “‘micro”’ viscosity of the solvent may be less than the ordinary 
viscosity. ‘The absorption band areas f «dw (Lewis and Kasha 1945) give about 
10-8 sec. for both solid and liquid anthracene, assuming the accuracy of the curves 
of Kortum and Finckh (1942). ‘The exciton movement in the crystal does not 
seem to cause any great lengthening of the period between absorption and radiation. 
‘The question then arises, what terminates the resonance transfer process in the 
absence of any added trapping agent? There are two possibilities : radiation may 
occur from molecules which are “‘ displaced” in the crystal, or at the surfaces of 
““mosaic structure’’, or the removal of excess vibrational energy from the excited 
molecule may stop transfer. Figure 2 shows potential energy curves for a diatomic 
molecule. 

Absorption is represented by the vertical line BK and fluorescence by GM. 
‘The time for passage from K to G is about 10-!sec. It might be assumed that 
exciton transfer continues while the energy available is represented by point K, 
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but that when a molecule falls to G it has insufficient energy for an energy exchange 
chain. . 

A further matter of interest is the alternative of radiation or quenching in these 
crystals. Quenching is undoubtedly favoured by polarizability effects (Bowen 
1947). The series of solids anthracene, naphthacene and pentacene have 

- fluorescent efficiencies at 20°c. of about 0-9, 0-002, and 0 respectively, probably 
due to increasing polarizability of the molecules. When embedded in an anthra- 
cene crystal naphthacene and pentacene molecules fluoresce with much higher 
efhiciencies. If two polarizable molecules approach closely, their potential-energy 
curves (see Figure 2) become shifted and less well defined, especially in the region 


POTENTIAL ENERGY 


DISTANCE 


Figure 2. 


marked X. Radiationless transition to the ground state is then possible, the 
excitation energy passing into vibrational energies of the two molecules. Reson- 
ance transfer of energy, however, occurs at larger distances, where polarizability 
has less effect, but where there is sufficient electronic interaction present. Both 
processes may be visualized as initiated by the transfer of an electron from an 
unexcited molecule to the vacated orbital of an excited one with transfer of the 
excited electron in the reverse direction. If the transfers are simultaneous the 
process gives rise to the exciton, but any check between them may lead to quench- 
ing. How these two effects depend on the relative orientations of molecules to 
each other with ensuing directional effects in the crystal is a matter about which we 
know little at present. 
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ABSTRACT. An adjustable Geiger-Miller counter and special electronic circuits have- 
been used to study the effect of several counter variables on the occurrence of multiple 
pulses. ‘The results show that the plateau characteristics combine the effects of an increase 
in efficiency and an increase in the probability of a multiple pulse as the operating potential 
is raised. It was found that the probability of occurrence of multiple pulses does not depend 
simply upon the charge per pulse, but that it depends also upon the anode diameter and the 
total pressure of the filling. 

Curves have been drawn showing the charge per pulse as a function of the operating 
potential. A change of slope occurs in these curves when the charge per pulse is approxi- 
mately equal to the charge on the active portion of the anode. 

The counter was used as an ionization chamber at operating potentials below the Geiger 
region. It was found that as the potential is raised above the threshold of the proportional 
region the ionization current varies with potential as log J oc V while at higher potentials the 
relation is adequately expressed by loglog Ja V. Near the Geiger region a sharp break 
occurs in the log log curve, which is taken to indicate that the discharge begins to spread at 
this stage. 


§1. INTRODUCTION 


N a previous paper (Chaudhri and Fenton 1948, subsequently referred to as I) 
it was shown that with a given cathode and filling the length and slope of the 
plateau depend strongly on the anode diameter. ‘The present paper is a 

report of work primarily intended to extend the investigations to include the effect 
of pressure on the plateau characteristics. 

Previous workers (see for example Korff and Present 1944, Curran and Rae 
1947) have attributed the slope of the plateau to spurious pulses and to an increase 
in the effective anode length as the operating potential is raised. It is convenient 
to divide the spurious pulses into two classes : (a) those pulses which occur as a 
result of previous pulses and are separated from them by time intervals of the order 
of several hundred microseconds; (b) pulses which are not the result of the entry 
of any ionizing particle into the effective volume of the counter. 

Pulses of type (a) may be distinguished by their non random distribution while 
pulses of type (b) can be distinguished if the rate at which they occur depends on 
the counter variables. In particular, to affect the plateau slope the rate must vary 
with operating potential. In view of the observation that the plateaux corrected 
for spurious pulses of type (a) are the same within experimental fluctuations over 
the range of pressures and anodes used, we believe that very few pulses of type (b) 
occur in our counter. 

To study the contribution of spurious pulses of type (a) to the plateau character- 


istics an interval discriminator circuit was used during the measurements. During 


the course of the work the charge per pulse was measured over the range of 
operating potentials for each anode and pressure. The counter was also used as an 
ionization chamber and the results provide interesting information about the 
threshold of the Geiger region. 
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For convenience we give below a list of the symbols and definitions used in 

the report. 

V, threshold of the Geiger region. 

M number of square waves per second produced by the multivibrator. 

C number of coincidences per second. 

R counting rate in counts/second excluding spurious pulses of type (a). 

S number of spurious pulses of type (a) per second. 

G total counting rate in counts per second. 

P_ measured probability that any pulse will produce as purious pulse of type (a).. 
P=SIG. 

ty dead time, i.e. the time which elapses after the occurrence of a pulse before. 
another Geiger discharge can occur. 

q _ charge per pulse. 

§2. APPARATUS AND EXPERIMENTAL METHODS 

The adjustable counter described previously (I) was used during the present 
observations. ‘The counter was set so that the central 5 cm. length of the anode 
was used, and the filling was 90° of spectroscopically pure argon plus 10% of 
alcohol. ‘The tungsten anodes used were of 0.05, 0-07, 0-1, 0-15, and 0:20 mm. 
diameter. ‘The gas pressure could be changed without pumping and refilling (a 
process which has proved to be not reproducible with our set-up) by adjusting the 
level of mercury in a 1-litre flask attached to the system. In view of the experi- 
ments reported by Korff (1944) no trouble was to be expected due to the small 
amounts of mercury vapour thus introduced into the counter. ‘Throughout the 
work to be described the pressure could be set at any value from 7:5 to 19cm. 
of mercury, the pressure being read from a U-tube manometer. Relatively 
low counting rates (of the order of 1,000 per minute) were used because 
earlier results had shown that high counting rates alter the background (1). 

To find the contribution of spurious counts to the plateau characteristics the 
circuits shown in Figure 1 were used with the triggered multivibrator set to give 
waves of 1,000 usec. duration, thus dividing the Geiger counts into two groups: 
group 1, recorded by scaler 1, counts occurring at any time greater than 1,000 psec. 
after the previous count, group 2 (scaler 2), those occurring less than 1,000 psec. 
after the previous count. 

A small correction was made for the number of random counts in group 2. 
Because of the dead time of the counter during which no pulses, random or other- 
wise, can occur, the effective time interval used for the calculation of the number of 
random pulses likely to occur is 1,000 usec. minus thedead time. The deadtimes 
were therefore measured, for each operating potential during the plateau measure- 
ments, by decreasing the length of the square wave until no coincidences occurred. 

It was shown that all spurious counts occur within 1,000 sec. of the previous 
pulse by measuring the rate at which spurious counts occurred for a series of values 
of the multivibrator interval at constant operating voltage. ‘The spurious rate 
increased rapidly with multivibrator interval at first but ultimately became more 
constant. ‘This measurement was made at points near the beginning, the middle, 
and the end of the plateaux for each anode using a pressure of 12cm. Hg. The 
time value at which the rate ceases to increase with interval was always found to be 
less than 1,000usec. A check was made at other pressures by increasing the 
interval from 1,000 to 1,360 psec. and as this gave no increase in Sit was concluded 
that all the spurious counts occur within the 1,000 psec. interval. 
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For the charge measurements a type 954 acorn valve was used in an electro- 
meter circuit attached to the anode of the counter as shown in Figure 1. In 
practice the opening of the key K was synchronized with the closing of the scaler 
counting switch, and the number of counts necessary to alter the potential of the 
grid of the electrometer valve by 0-5 volts was registered. Hence the average 
charge per pulse can be calculated, using the measured capacity of the system. 


Scaler | 


Differentiator 


Triggered 
Multtibrator 
20 to 360,usec: 


Cathode 
1 Follower 
Amplifier 


Amplifier 


Arrangement for plateau measurements 


Scaler 2 


Electrometer 


Cathode 
Follower 
Jk 
‘ Amplifier 
OE.H.T. + 


Arrangement. for pulse size measurements 
Figure 1. 


When the counter was used as an ionization chamber the electrometer was 
employed to indicate the collection of a predetermined charge in a measured time. 
By altering the capacity C it was possible to measure over a large range of ionization 
currents without having to measure either very short or very long time intervals. 
The currents measured ranged from about 2 x 10-15 to 10-§ amp. when the gamma- 
ray source was placed close to the counter. An additional factor of about 100 was 
obtained in two stages by moving the source away from the counter, so that in 
effect the measurements covered the current range from about 2 x 10-18 to 10-8 amp. 


Sis, JASON eh SS) 
3.1. Plateau Characteristics 


The curves obtained have the general form shown in Figure 2. The corrected 
plateaux corresponding to the various pressures and anode diameters give the 
shape of the plateaux that would be obtained if there were no multiple pulses 
Graphs of these lie close together when plotted with the same counting rate and 
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jj overvoltage coordinates. As there are no observable trends ‘of variation with 
| pressure or anode diameter it appears that the shape of the corrected plateau is, 
from the practical viewpoint, independent of these variables over the range covered 
| by the present experiments. Eventually, as the potential is raised, the corrected 
i plateau breaks off due to sprays of multiple pulses lasting more than 1,000 psec. 
The curves of P plotted against operating potential vary considerably with 


Co 
Cs) : 
<= 5 
Seah ] 
= 
=) 
= 
~ 3-0F 
mo) 
Cc 
3 
22:5 = 
i= 
a 
ES ) 
20+ 
5 40-75 S 
oO 
2 ISP | & 
= 3 
(2 lls 
2 10r 5 
8 s 
n dn. wo 
—— 4 ari 1 3 OSPF ite = 
‘1400 1500 1600 1700 1800 2 Se 
Anode Potential (volts) ea 0 fl : : ; <. 
! : 0 100 200 300 400 500 600 700 800 900 
igure 2. Typical plateau analysis. Multivibrator Interval 7’ (z sec.) 


Figure 4. 
Counter data: Anode diameter 0-07 mm. 
Total pressure 12 cm. 
Operating Potential: A and B 1,440 v. 
Cand D 1,240 v. 
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3. Variation of P with overvoltage tor Operating Potential (volts) 


' various anode diameters (d) and total 
_ pressures (7). Figure 5, 


pressure and with anode diameter (Figure 3). For purposes of comparison it is 
convenient to plot, as a function of the potential, the probability that a multiple 
count occurs, instead of the spurious counting rate. ‘The results show that P 
increases rapidly at overvoltages corresponding to the upper ends of the Geiger 
plateaux; the overvoltage at which this rapid increase begins increases as the 


Sam 
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anode diameter is reduced from 0:2 mm., passes through a maximum at an anode 
diameter of about 0-1mm. and then decreases. The slope of the straight 
portion of P curves increases the anode diameter is reduced. 


3.2. Spurious Count Distribution 


The curves obtained by plotting the spurious counting rate against multivibrator 
interval show that the spurious counts begin to occur immediately after the dead 
time (Figure 4, curves Aand C). As measured above, f, is the time after a pulse 
during which the counter is unable to give a pulse as great as the minimum detectable 
size. By observing the pulses on a triggered cathode-ray oscillograph it was 
shown that this minimum size is reached after about 1-1 of the true dead time and 
that pulses smaller than this minimum detectable size also occur more frequently 
than the random rate. 

The distribution in time of the spurious counts was obtained by taking a 
histogram of the curves of S plotted against multivibrator interval. Except at 
high overvoltages the maxima of these histograms occur at a time approximately 
equal to that required for the space charge to cross the counter, assuming an ionic 
_ mobility intermediate between those given in the International Critical Tables 
for argon (k=1:66cm/sec. per v/cm.) and alcohol (k=0-48 cm/sec. per v/cm.). 
Figure 4, curve B, is for a high overvoltage. At overvoltages near the middle of 
the plateau and below the histograms are nearly symmetrical (Figure 4, curve D). 


3.3. Pulse Size Measurements 


Figure 5 shows the average charge per pulse g plotted against the operating 
voltage for a series of pressures using the 0-1-mm. diameter anode. The line 
q=Qzy gives the charge on the active length (5cm.) of the counter wire. With 
anodes of diameter 0-05, 0-07 and 0-1 mm. and at pressures ranging from 8 cm. to 
18cm. Hg a sharp break occurs in the curves at about g = Q) while the ratio of the 
slopes of the lines on either side of the break ranges from 0-4 to 0:67. For larger 
anodes (0-15 and 0:2 mm. diameter) and at the higher pressures the break occurs 
for g<Qy and the slope ratio increases. 

In Figure 6 the charge obtained at a total pressure of 10cm. Hg is plotted against 
overvoltage for four different anode wires. ‘The curves show that the charge at a 
given overvoltage increases slowly with increasing anode diameter. 

The break voltage also increases with increasing wire diameter. ‘The starting 
voltage, from which these overvoltages are obtained, is defined as the voltage at 
which the extrapolated q~V curve shows g=0. Figure 7 shows that the break. 
voltage increases linearly with pressure for all anodes where the break is clearly 
defined. 


3.4. Measurement of Ionization Current 


‘The counter was used as an ionization chamber, connected as shown in Figure 1. 
As the potential is raised above the proportional starting voltage V, the 
logarithm of the ionization current J is proportional to the potential at first ; 
then at higher potentials log log J becomes proportional to potential as shown in 
Figure 8, curves Aand B. Sharp breaks occur in the log log curves and when 
these are matched, as shown in the figure, the curves for different gas pressures 
converge as the potential is raised. Figure 8, curve C shows J plotted against 


potential for the region above the break for 16cm. pressure. This curve may be 
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compared with the charge per pulse curves of Fi 

gure 5. Figure 8, curve C 
shows that the value of V, as defined in §3.3 is about 50v. above the potential 
: igs by the break in the log log curve. Figure 9 shows that the value of 
| log log J at the break in the log log curves depends linearly on the pressure. 
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3.5. The Dependence of P on q 


The probability P of the occurrence of spurious pulses and the charge per pulse 
qg were obtained as functions of the anode length, and from the resulting straight 


line relationships the relation between P and q was obtained. Figure 10 shows the ~ 


linear dependence of P on qg for a 0-1mm. anode, 10% alcohol, and 12cm. Hg 
total pressure, at a fixed operating potential of 1,400v. This method of varying 
q has the advantage that as the operating potential is held constant the discharge 
mechanism does not alter, the charge increasing with anode length simply because 
the discharge has further to spread. ‘The possibility that the variation of counting 
rate with anode length may have some effect on the value of P was removed by 
using a fixed anode length and altering the position of the gamma-ray source to 
cover the range of counting rates obtained during the experiment in which 
the anode length was varied. P was found to be independent of R over this range. 

In order to find the dependence of P on the anode diameter and the pressure of 


the filling a fixed value of g was chosen and the corresponding values of P were | 


taken from the Pvs(V—V,) curves, examples of which are shown in Figure 3. 
The results are shown in Figure 11 where the values of P corresponding to 
q=8 x 10-?° coulombs are plotted against pressure (curve A) and anode diameter 
(curve B).. This value of g was chosen because, with the exception of the value at 
8 cm. total pressure the corresponding values of P lie on the early portions of the 
Pvs (V —V,) curves, in the region where the majority of the multiple pulses contain 
only two individual pulses. The rapid increase of P corresponding to the upper 
limit of the plateau is associated with multiple pulses of multiplicity greater than 
two. 


§4. DISCUSSION 


The corrected plateaux of which an example is given in Figure 2 do not flatten 
to any extent until the operating potential is appreciably above the threshold. 
Further, as can be seen with a cathode-ray oscillograph, the pulses are of approxi- 
mately equal size before the apparatus begins to register counts (minimum 
detectable pulse size=1-2v.). It has also been shown by varying the amplifier 
gain that the increase in the corrected counting rate with operating potential is not 
due to small pulses being brought within range of the apparatus. Part of the 
increase may be due to an increase in the effective length as mentioned in § 1, but 
because of the method of shielding the anode used in our counter we do not think 
that end effects can account entirely for the observed increase. On the other hand, 
plateaux of this form are predicted by Narwijn (1943) and Wilkinson * who obtain 
the probability that the discharge begins to spread as a function of overvoltage. 

For one electron set free in the counter by the ionizing agency this probability 
is rather low (0-6) at an overvoltage of about 50v. (i.e. when counts are first 
registered) and increases with overvoltage to about unity at 200 v. above the thres- 
hold, corresponding to a fairly large slope in the corrected plateau. At a fixed 
overvoltage the probability of spreading increases as the specific ionization in- 
creases, and at 50 v. overvoltage is practically unity when the number of electrons 
set free by the ionizing particle exceeds five. A gamma-ray photon can only be 
detected if it releases an electron which then enters the effective volume of the 
counter. In our case the gamma rays from radium are filtered by 3 cm. of lead so 
that most of the primary electrons will have energies of 1 Mev. or over. Depending 


* We are indebted to Dr. D. H. Wilkinson, Cavendish Laboratory, Cambridge, for this inform- 
ation prior to its publication. 


a ee oa 


4 


Experiments with an Adjustable Geiger- Miller Counter 39 


on its path through the counter, from one to about eight secondaries will be pro- 
duced by an electron of this energy, so that we should not expect the plateau 
to flatten until appreciable overvoltages are used. From the practical standpoint 
the corrected plateau curves show that in a counter such as ours the overvoltage. 
used should be of the order of 200v. This would apply also when the cathode is 
treated to reduce spurious pulses to a minimum, or when a paralysis circuit is used 
to prevent their being counted. 

It was reported earlier (I) that there is an optimum anode diameter correspond-: 
ing to a given filling and cathode for which the best Geiger plateau is obtained. 
This observation is explained in terms of the way in which the spurious counting 
rate curves vary with anode diameter and the way in which they combine with the 
corrected plateaux to give the Geiger plateaux. It should also be pointed out that 
the sum of the corrected plateau and the spurious counting rate curve does not 
always yield a curve with a really straight portion. This explains the deviation of 
points sometimes observed from the straight lines usually drawn for the Geiger 
plateaux, even when a large number of counts are taken in order to reduce the 
standard deviation. 

The sharp breaks in curves A and B of Figure 8 are interpreted as marking the 
stage at which the discharge begins to spread along the anode. This view is 
supported by observations with a cathode-ray oscillograph in the neighbourhood 
of the break. A few volts above the potential corresponding to the break there is a 
distribution of pulse sizes. The maximum pulse size observed at a given 
operating potential is halved by halving the effective length of the counter, showing 
that at quite low overvoltages the discharge can spread to the full length (10 cm.).. 
As the potential is raised the number of pulses reaching the corresponding maxi- 
mum size increases until at about 20 volts overvoltage all the pulses appear to be of 
the same size. The starting potential for the Geiger region is generally defined 
(Korff 1944) as the potential at which all the pulses are of equal size. While this 
definition is usually adequate for the practical use of counters it simply indicates. 
the potential at which the discharge spreads to the full length of the counter, 
and because of the statistical nature of the spreading mechanism (Alder et al. 1947,, 
Wilkinson, unpublished) its value depends on the counter length. However,. 
when the discharge mechanism is being investigated the potential at which 
spreading begins is often of importance and this may be determined quite simply by 
the method described in this paper. 

The observed dependence of the probability of occurrence of spurious counts on 
the charge per pulse is consistent with the view expressed by Korff (1944) that. 
multiple pulses are due to statistical breakdown in the quenching mechanism 
which occurs when the positive ions reach the cathode. However, P does not 
depend only on the number of ions reaching the cathode but, as shownin Figure 11, 
depends also on the anode diameter and the pressure. These results suggest that 
the ions which arrive at the cathode are not all of the same type, and that the 
production of the ions which are more likely to produce spurious pulses is favoured 
by the use of high pressures and small anode wires. The effect may be due to the 
increase in the field near the anode when smaller anodes or higher pressures are 
used. Calculation of the fields at the anode surfaces using the operating potentials. 
corresponding to the points given in Figure 11 yields values of the field strength 
ranging from 37,000 v/cm. for the 0:2 mm. anode to 110,000 v/cm. for the 0:05 mm. 
anode. For the 0:1 mm. anode it ranges from 49,000 v/cm. at 8cm. pressure to 
64,000 v/cm. at 18cm. pressure. 
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The spurious count distribution curves (e.g. Figure 4) together with the 
-observations with a triggered cathode-ray oscillograph show that some spurious 
pulses occur as soon as the end of the dead time is reached and that the pulses are 
distributed over a fairly wide range of times. These curves also provide evidence 
that the ions formed in the discharge are not all of the same type. ‘The distri- 
butions obtained at low voltages are symmetrical, as would be expected if the space 
charge is regarded as a group of ions moving with a Maxwellian velocity distri- 
bution upon which one drift velocity (due to the electric field) is imposed. This 
applies in the outer part of the counter where the field varies slowly so that 
focusing is negligible and ions of the one type have approximately constant {rift 
velocity. At high voltages the distribution curves are characterized by a sharp 
rise to the maximum followed by a very gradual fall. Such a distribution would 
be given if the space charge contained multiply, as well as singly ionized atoms; 
the higher drift velocity imposed on the multiply ionized group would give the 
early maximum and the subsequent arrival of the singly ionized group would 
explain the very gradual fall. This view is supported by the fact that the asym- 
metry sets in at points lower on the plateau for the smaller anodes, as here the 
field strength at the wire is greater. The continuous electron transfer between 
singly ionized argon atoms and alcohol molecules which, according to Korff (1944), 
explains the quenching action, does not affect the argument. Such a transfer 
affects the position of the maximum because of the difference in the mobilities of 
the ions, but it could not produce the observed asymmetry as the exchange is 
complete before the ion sheath reaches the cathode. 
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§ ABSTRACT. The Smoluchowski equation for the distribution of the relative coordinates 
of a representative pair of molecules in a liquid is formulated in such a way that the effect 
of non-uniform velocity or temperature is taken into account. The equation is reduced 
to an ordinary linear differential equation of the second order. One of the two constants 
of integration is readily determined whereas the determination of the other constant may 
lead to inconsistencies. The molecules are assumed to interact with the average potential 

that is appropriate to thermal equilibrium. If it is assumed that this potential becomes 
infinite of the first order at distance zero then the problem of viscous flow and the problem 
of thermal conduction are uniquely soluble. Approximate evaluation of the transport 
coefficients shows the right kind of temperature effect. “lhe result is in accordance with 
the previous treatment by means of the cell model. 


§1. INTRODUCTION 


N a previous paper (1947) referred to as I the author analysed the transport 

phenomena in liquids by means of extremely simplified assumptions on the 

liquid state. The “cell model’? was employed and the Smoluchowski 
equation was used for determining the distribution in coordinates of one repre- 
sentative molecule. The influence of the rate of shear or gradient of temperature 
-on the distribution was taken into account only on the surface but not in the 
interior of the cell. The calculated distribution functions were infinite at the 
centre of the cell but yielded finite values for the flow of momentum or energy. 
_ The marked difference of the temperature effect on viscosity and thermal con- 
-ductivity was obtained without any additional assumption and in spite of the 
complete similarity of mathematical operations which lead to the one and 
to the other transport coefficient. 

The author was unaware of a paper by Auluck and Kothari (1944) in which 
formulae similar to the author’s are derived from the well-known “ hole” model 
as given by Fiirth (1941); in addition, numerical values of the right order of 
magnitude are obtained. In their paper the mechanism and the rate of transfer 
is explicitly assumed to be markedly different for momentum and energy. ‘The 
author’s result shows at least that assumptions of this kind are unnecessary. 

In the meantime Born and Green (1947) published the extension of their 
general kinetic theory of liquids to transport phenomena and made a preliminary 
calculation of the viscosity; in their result various approximations are involved 
and the equilibrium distribution of a pair of molecules in a liquid is replaced 
by the distribution of a pair of isolated molecules. Since thermal conductivity 
still awaits detailed treatment, and the introduction of the liquid state distribution 
‘may affect the result appreciably, it cannot yet be seen to which conclusions the 
rigorous theory, as developed by these authors, will lead. 

The method as employed in paper I is too crude for the actual calculation of 
the transport coefficients. In addition, its foundations, that is, the cell model 
.and the application of the Smoluchowski equation, are open to criticism. 
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The cell model is hardly suitable for a quantitative theory since it is impossible 
to divide the liquid into cells of fixed size and shape without introducing arbitrary 
and sometimes inconsistent assumptions. The partial success of this model 
is due to the fact that fundamentally it is an approach to the specification of the 
liquid state in terms of the distribution of relative coordinates of two representative 
molecules. : 

The Smoluchowski equation is at first sight a legitimate approximation to 
the equation which has been shown by Kirkwood (1946) to be adequate and it is 
preferable from the point of view of mathematical simplicity. Its application 
in paper I resulted, however, in distributions which are infinite at one point. 

In the present paper an attempt is made to clear the ground for the calculation 
of transport coefficients by employing the pair distribution instead of the cell 
model and by admitting only solutions of the Smoluchowski equation which, 
when interpreted as probability distributions, are in agreement with the physical 
conditions. 


$2. THE MODE EY OR Eh s 1@ Gil Di Sue AGE: 


In this paper the liquid state is described in terms of the relative coordinates 
of two representative molecules. The probability distribution of the relative | 
coordinates of the two molecules is, in the case of monatomic liquids, a probability 
distribution of the distance ry, The assumption is made that the two molecules 
interact with a force equal to the negative gradient of an average potential O(7) 
which itself is derived from the equilibrium distribution. 

The probability density in equilibrium is zero at small distances, has one 
strong and possibly several weak maxima and tends to a constant value at large 
distance. It can be derived from experimental x-ray interference patterns or 
calculated from the intermolecular forces according to the theory of Kirkwood 
(1935) or Green (1947). 

The average potential ® which is essential for the model under consideration 
is defined by writing the equilibrium distribution g(r) as 


a(ry=(1/VoexpC— O/R EW 0 ek eee (1) 


® depends according to this definition explicitly on the temperature, but little 
is known of the amount of this dependence. In equilibrium only the minima 
of potential are densely occupied; in non-equilibrium distributions all ranges 
of the potential curve may be of importance. For large distances ® is zero. 
Where the potential is markedly positive it is not likely to depend very much on 
the temperature. Near the main minimum at about one molecular diameter 
distance, the potential varies with temperature so that the minimum becomes 
flatter when temperature decreases; this variation is an expression of the fact 
that the molecules have not the tendency to crowd together at the distance of 
the minimum. At small distances the potential is likely to approach the repulsion 
potential of a pair of isolated molecules. If overlapping of the electronic clouds 
is almost complete—the corresponding small distances are practically never 
reached—the potential becomes infinite of the first order owing to the Coulomb: 
repulsion of the nuclei. 

For the model. under consideration the Smoluchowski equation has to be 
modified in order to account for the non-uniform velocity or temperature inde- 
pendently of the boundary conditions. 
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§3. THE DIFFERENTIAL EQUATIONS 


| If temperature and velocity are uniform, the distribution, g, of relative co-- 
| ordinates is determined by the Smoluchowski equation 


Og/ot =div[grad(kTg/Bm)+(ggrad@/Bm)], — ...... (2) 


where f is the “friction constant” having the dimension time, — grad O/Bm 
is the mean displacement, k7'/8m is the mean square displacement per time unit ;. 
the expression in the bracket is equal to the negative flow vector j. In thermal. 
equilibrium g =g,, j=0. 
. For a steady non-uniform state let g=g,(1+w). For the specification of 
viscous flow let the velocity ¢ depend in an arbitrary way upon the coordinates 
| subject to the condition dive=0; ce is added to the mean displacement. 
| Neglecting terms of the magnitude c.w, the flow is then given by 


toe (RL pm) oradm |e ee (3 a) 
and w is determined by 
V?w —(grad Y, gradw)= —(Bm/kT\(c, grad VY), (4a) 
where ‘’=®/kT subject to the boundary condition 
radii (e (7 1ce)" | nu cc, (5 a) 


so that in infinity the flow is equal to gy. ¢. 

For specifying thermal conduction let T be an arbitrary function of the- 
coordinates subject to the condition V?77=0. Neglecting terms of the magnitude 
w grad T, grad? T etc., the flow is given by 


j= —(Rk/Pm)g,[T,) gradwt(1+'V)gradT], _...... (3d) 
where T is the temperature at the origin and w is determined by 
V2w —(grad V, gradw)=(¥/7,)(grad VY, grad T), ...... (4) 
subject to the boundary condition 
Lygracwe: Cratt Oe esate (5 b) 


so that the flow is zero at infinity. 

For the purpose of calculating viscosity let the flow be laminar with a rate of 
‘shear equal to 4a and the velocity be distributed according to c,=c,=0; c,=ax.. 
Then the right hand side of (4a) is —«¥’rsin?@ sing cosf; «=aBm/kT. Let 
w have the form w=u(r)sin?6 sing cos¢. The function w changes sign at. 
f=0, 7/2, 7, 37/2. If one molecule is placed arbitrarily at the centre the pro- 
bability that the other will be in one pair of quadrants, say 1 and 3, is greater’ 
than that it will be in the other quadrants 2 and 4. This means that in the 
flowing liquid the distribution of molecules is anisotropic. The function w is- 
determined by 


u” +[(2/r) —V" Ju’ —(6u/r?) = —al'7, ne ees (6a) 
with the boundary condition 
dujdy =O UiraO O(r= GO) Nileaee (7 a) 


For the purpose of calculating thermal conductivity let the temperature be 
distributed according to T=7,(1+6z) so that the gradient 1s constant: 
grad, T=grad,T=0; grad, T=67). Then the right-hand side of (46) is. 
bY'L’ cos@. Let w have the form w=v(r)cos@, then v is determined by 

vo” +[(2/r) —W"]o’ —(20/r?)=BY eee (6) 
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with the boundary condition 
jidr=a bs er bs (FSH). 1 ee (7b) 
The calculation of the transport coefficients is reduced to the integration of 
the radial differential equations (6 a), (6). Which solutions are admissible must 
be decided by considering how the probability for finding the relative coordinates 
in a finite volume deviates from the equilibrium distribution. This deviation of 
probability is proportional to 


Ja= Lour® dr and J,=| govr’ dr 
am 0 


respectively. These integrals must be finite for every finite value of the upper | 
limits. 

It is known from experiments that laminar flow does not cause any appreciable 
anisotropy in liquids composed of simple molecules. The integral J, to which 
the probability of preferential distribution in an odd quadrant is proportional, 
‘gives a measure of anisotropy, and it follows from what has been said that even 
if the upper limit becomes infinite the integral J, must remain finite. No 
restriction of this kind is expected to hold for the integral J,. 


§4. SOLUTIONS NEAR THE SINGULAR POINTS 


Integration of equations (6a), (65) near their singular points at the origin 
cand at infinity provides the means for selecting solutions of physical significance. 

Y has a pole of the first order at the origin; it will alternatively also be assumed 
ithat ‘V(0) is finite. ‘Y’ =0 if 7 is sufficiently large. 

The two arbitrary constants of integration are denoted by A,, A, and B,, By. 
‘Other constants appearing in the solutions are denoted by F and H with suitable 
-sufhixes. 

The homogeneous equation corresponding to (6a) is 


uo (Zin) Ee = (Gur ye i ae (8 a) 
If the term — ‘fw’ is omitted an equation results which has the solutions «=7? 
and u=1/r°. 
At large values of r equation (8 a) has accordingly two solutions, one being 
finite at the origin and of the form 


Wy eal aka TOUS) el a cee (9a) 
the other vanishing at infinity and of the form 
yore OF at gael ak ges orn sh (10 a) 


Assuming at first that ‘’(0) is finite, it is seen that at r=0 the expansion of u, 
‘commences with the second power of r; it can be written as 


ae ceca os 4 (|) ees a eh (11a) 
‘The other solution has near the origin the form 
Ug= Fejr+.... Feed Wr) 
The general solution of equation (6 a) can be written as 
usAid, PAgs tis = 3 ee ae (13 a) 
where Us—tor, voaees (14a) 


is a particular solution. 
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If the potential has a pole at the origin so that for r <7, it is equal to Y =C/r, 
this does not affect w, and the expressions (9a), (10a). In the range 0<7r<1, 
a the appropriate expressions for u,, uw, are found to be 

Dies u,* =1 +4 (6r/C) + (12r2/C2); i 
ty = Fity* + Feig* 5 uy* =[1—(6r/C) + (1272/C2)Jexp(C/r) | 
| Equation (6d) is discussed in a similar way. The corresponding homogeneous 

- equation 


| Ut [(2r/v—¥"]o’ —Qo/r*N=0 aes (8b) 


il y has, for large values of r, solutions of the form 
= Orca e bye aye GY) gale AN, (emi) Met aead ON lia dN he rg ops (9b) 
V_=1/r? + 0(1r-%). Se esaees (EOD): 
_ If (0) is finite the expansion of v, at r=0 can be written 
Detisca’s 1 AO) Roleemeene 1 eee (11): 


The other solution has near the origin the form v,=H3/7? + . . 
The general solution of (6 5) can be written as 


Pi iseads Uigats IsgU peter as te. aw Al eis atirees (13 4). 
where the particular solution v, is given by 
v3(r) = — b| POY Oexp{-— FOF [ear —vu(rolO dé, (16) 
and has at large values of r the form 
v, = bH or + bH,/r? + 0(7-). weeeee (146) 


If the potential has a pole at the origin the appropriate expressions for 2, V2, 
in the region 0<r<7, are found to be 


Peay Cine Gl, \ Bane: (15d) 
vg = Hyv,* + Hgv2*;  v2* =[1—(2r/C)]exp(C/r) | 
| and v,* =b{[1 + (27r/C)][2(79 —7) + C In (r/79)] 
: +[1=(2r/C)]f2r—Cexp(Cir)| exp(—2)E ae} 


is substituted for v,. At the origin the function v,* and its derivative 
doar |(Ci7*) + (2/7) —(2/C)] exp (Cir) =. eee (17) 


are exponentially infinite and v, is infinite of the order Inv. 


§5. SELECTION OF SIGNIFICANT SOLUTIONS 

In determining the constants of integration it is at first assumed that ‘Y’(0) 
is finite. 

In the problem of viscous flow the first condition (7a) is applied to (13 a); 
considering (9 a), (10 a), (14a) an equation determining A, is derived: 

PASH eth )a0 tA iin eh ti) SPV It ie ioners ie: (18 a) 

The corresponding choice of A, complies with the second condition (7 a). The 
constant A, has to be chosen so that the probability J, is finite for infinite 7. 
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If A, is assumed to vanish, J, becomes infinite owing to the second term in (9 a). 
If A,~0, the integral determining J, is divergent. No solution of (6a) is 
admissible. 

In the problem of thermal conduction the first condition (76) is applied to 
(136); considering (96), (106), (146) an equation determining B, is derived: 


DUTT ACh iio. inal ae eee (18) 


“The corresponding choice of B, complies with the second condition (7b). ‘The 
integral determining J, converges independently of any choice of the constant 
B, although the integrand is infinite at the crigin if B,~40. The assumption 
B,=0 is most plausible but not necessary; the distribution and consequently 
the thermal conductivity remain undetermined. 

These difficulties are overcome by letting the potential have a pole at the 
origin. In the problem of viscous flow A, is determined so that r?u vanishes 
-at infinity : 
Ay Bs Ag Qn | man ae Sete (19 a) 


Whereas u is now exponentially infinite at the origin, the integrand of //, is finite 
and the integral is finite for infinite distance. In the problem of thermal con- 
-duction the constant B, is not determined by considering J,. The integrand 
of (256) determining the flow of energy contains terms which are proportional 
to (rt —‘V)gyv'r?. At the origin the expression in the bracket is proportional 
to 1/r and according to (155) and (17) wv’ is proportional to 1/r? if there is any 
‘contribution of v,* to v. In order to prevent the flow of energy becoming 
infinite it is therefore necessary to let 


Byx0, “SP (2.5 o6 Oe Sees (19d) 


The necessity of specifying the potential near the origin in detail is due to 
the limit of validity of the Smoluchowski equation. Near the origin the calculated 
functions cannot be expected to be good approximations; in the problem of 
viscous flow the correct probability density is likely to be much smaller than 
gou. The contribution of the region near the origin to the flow of momentum 
or energy is, however, small since the corresponding integrands vanish. 


§6. APPROXIMATE INTEGRATION 


Approximate expressions for the distributions are found by solving differential 
equations which are distinct from the correct equations by relatively small terms. 
If ‘V(0) is finite and y is a constant let 


flr)= 3? yr [exp Meds, 
‘This function solves the ejuation 
u" +[(2/r) —P"]u’ —(6u/r?) — (3, —( —y)r? (Pu /rf,) =0, ...... (21a) 
which differs from (8 a) by the last term only. ‘Transformation of (8 a) by means 
of 
Url exp the 7 Pee aah en ee (22 a) 
‘gives Uy + (CE 7) Pa = EP" (6/7?) 0 ee ee ee (23 a) 
Substituting (22a) in (21a) with y=0 and y=1 yields equations differing from 
(23 a) by 2\'u9/r and —3"’u)/r respectively. ‘They are small compared with 


a 
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4 | the second term of (234), particularly if |'¥’|>1/r or ¥’=0. f, is therefore 
“an approximation to w, and the other solution of (21 a) 


fii, PO | [ea Oltexp FOds pea (24a) 


) is an approximation to uy. 
A value between 0 and 1 is then assigned to the constant y in order to obtain 


pie.) 
y 


a the best possible approximation at small values of r. The first two terms in the 


® expansion of f, are accordingly fitted to the corresponding terms in (11a). It 
@ is found that y= § exp(—'‘¥(0)). 

} In order to account for the pole of ‘ at the origin, f, and /, are calculated 
as appropriate to a modified potential curve which follows the tangent of the real 
9 curve at 7) and has the value 2C/r, at r=0. At r=r, the functions f, and u,* 


ip are joined and similarly the functions f, and u,* : 


SOT <7, #4 —u,*, u,= Fyu,* +P u*; -forr 1, uy = Ff, + Fofoy 9tts=So- 
The four constants are determined by four equations expressing continuity of 
the functions and their derivatives at 79. 

Approximations to the solutions of (85) are found similarly by means of the 
equation 


v’ +[(2/r) — PJ’ — (2v/r?) —[2h, —(1 —y)r] (8 e/rhy) =0, ...... (21 5) 

the solutions of which are 
h,(r) =4r+ yr| fexph(c) ete jaw ae Wiese (205) 
hy=ghy  g(r)=—| [eh Ob expV@Ode es (248) 


with y=} exp(—‘Y(0)), and the pole is accounted for in the same way as above. 
‘The function v, is equal to v,* for r<7, and is determined by (16) for r>79 with 
7, instead of 0 as lower limit of integration. 

The present approximations, even without corrections for the pole, are sufh- 
cient for estimating the temperature effect. When better approximations are 
required they can be obtained by substituting the first approximation for the 
function in the last term of (21 a) and (216), adding the modified term to the right 
hand side and integrating by means of standard formulae. 


S7: GHEE LEMPERATURE ERFECT 


The coefficient of viscosity 7 and the coefficient of thermal conductivity A 
are determined by expressions corresponding to those given in §4 of paper I: 


n= —(N°RT/a V)| go’ Sin > oisin7 Zod? du dds san (25 a) 


A= —(N?K2T/2BmV) { (FP —P)[(o'/b) + ¥ + 1] cos? 8 


+ (drb’ — P)[(v/rb) + ¥ +1] sin? O}gor? sin @drdbdp. —...... (25 b) 
In (256) Ty is replaced by T. 

Viscosity is known to have a factor exp(U/RT) where U is an energy of the 
magnitude of the potential barrier outside the main minimum and exp (U/RT) >1. 
It is accordingly of interest to find out which solutions of (6a) and (6 6) are of 
a magnitude similar to that of the exponential. 
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The integrand in (20a) and (206) determining f, and h, is small at small 
values of r but increases exponentially outside the main minimum. Accordingly — 
at large values of r the integrals are of the same magnitude as exp(U/RT). A | 
lower limit for f, is derived from (20a): : 


Nae ith 7 | fide =a Lyi ewe (26) 


and a similar expression applies to hy. | 

The functions f, and h, are, according to (24a), (246) extremely small at 
distances outside the barrier and reach the magnitude of unity only near the 
barrier. ‘Their possible increase at the minimum is limited by the lower limit 
of f, and h, according to (26), so that their value is of a magnitude smaller than 
that of the exponential. Only at small values of 7 will the functions f, and h, 
increase exponentially. 

The integrands in expression (16) determining the function v, have at no 
point the magnitude of exp(U/RT) so that the integrals do not reach this magni- 
tude; vs itself is of the magnitude of the exponential at large values of r owing _ 
to the factor h,(7). 

In the integrands of (25a) and (255) the functions u,, v, etc., are multiplied 
into gj). The high values of uw, and v, near the origin and of the functions 
and v, near the barrier are offset by the negative exponential in gj. At large dis- 
tances the high values of the functions wu, and v, are ineffective because the 
integrand vanishes with ‘”. 

The transport coefficients can have accordingly an exponential temperature 
factor only if a constant of integration is proportional to a factor of this kind. _ 
These constants are determined by the constants in (9a), (96). According 
to (20a), (206) and considering the magnitude of the integrals in these 
expressions, only the constants F, and H, have the magnitude exp(U/RT). The 
constants F, and H, and, according to (146) and (16), the constant H, are of 
smaller magnitude which by way of comparison may be said to be of magnitude 
unity. Determination of A, according to the boundary condition (18a) and of 
B, according to conditions (18) shows that these constants have no exponential 
factor. B, vanishes according to (196). Only the constant A, as determined 
by (19 a@) has an exponential factor which is peculiar to the viscosity of liquids. 
The characteristic difference of the temperature effect on the two transport 
coefficients is seen to be represented by the constants of integration. 


§8. CONCLUSION 


The effect of temperature on the transport coefficients as obtained in this 
paper is of a similar kind to that derived in paper I. The previous interpretation 
of the temperature effect is confirmed. In viscous flow the molecules are forced 
to occupy frequently positions of high potential energy whereas no such necessity 
is imposed on the molecules in the case of non-uniform temperature. In the 
first case the distortion of the equilibrium distribution is of higher magnitude 
than in the second. 

The methods applied in this paper and in paper I may be regarded as different 
approximations to the exact theory. ‘The present method no longer depends 
on the crude parabolic potential previously employed and can be adapted to the 
correct potential curve. 
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| ABSTRACT. A theory of yielding and strain ageing of iron, based on the segregation of 
carbon atoms to form atmospheres round dislocations, is developed. The form of an 
atmosphere is discussed and the force needed to release a dislocation from its atmosphere is 
roughly estimated and found to be reasonable. The dependence on temperature of the 
|— yield point is explained on the assumption that thermal fluctuations enable small dislocation 
loops to break away ; these loops subsequently extend and cause yielding to develop. 
catastrophically by helping other dislocations to break away. The predicted form of the 
relation between yield point and temperature agrees closely with experiment. 

Strain ageing is interpreted as the migration of carbon atoms to free dislocations. \ The 
| rate of ageing depends upon the concentration of carbon in solution and the estimated 
| initial rate agrees with experiment on the assumption that about 0-003% by weight of 

carbon is present in solution. 


$1 INTRODUCTION 


T was shown recently (Cottrell 1948) that the relief of stresses round a disloca- 
tion by foreigmatoms in solid solution can cause an equilibrium “‘ atmosphere ” 
to form, in which large solute atoms are gathered in the dilated part of the 

dislocation field and small ones in the compressed part. An atom in substitutional 
solution relieves hydrostatic stresses and its interaction energy with a positive 
edge dislocation has been given as 


ey ee a (1) 
3 i Wea ae 2 
where 7 and « are its coordinates relative to the dislocation, « being measured 
from the slip direction, G and v are the rigidity modulus and Poisson’s ratio 
respectively, is the slip distance in the dislocation, and 7, and 1,(1 +e) are the 
respective atomic radii of solvent and solute. This formula fails at the origin 
because the elastic continuum theory is not valid in this region. To a first 
approximation, the dilatation is zero at the centre of a dislocation andit is reasonable 
to assume that V =0 whenr<7rp, where 7)~A, and is otherwise given by equation (1). 
In a dilute atmosphere the concentration of solute will be 
MG Ne ex Pia (Tr a)iR Tj ene 6 geet (2) 
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where my, is the average concentration expressed as atoms per unit volume, but 
in a dense atmosphere substantial relaxation of the hydrostatic stresses is possible 
and an upper (saturation) limit of concentration is set by the condition that the 


' local dilatations produced by the solute atoms and the field of the dislocation are — 


everywhere equal and of the same sign. 

Dislocations surrounded by atmospheres can produce plastic flow in two ways. 
If the applied force is small, the dislocations cannot escape from their atmospheres 
and the solute atoms must migrate with the dislocations. By applying a sufh- 


ciently large force, however, the dislocations can be torn from their atmospheres, __ 


in which case they then become highly mobile and able to produce rapid flow 
under smaller forces; a sharp upper yield point, followed by flow at a lower yield 
point, occurs. A specimen which is unloaded in this overstrained condition 
contains free dislocations and, on immediate reloading, shows no yield point, 


but if it is rested for a sufficient time at not too low a temperature before reloading, 
solute atoms migrate to the dislocations to form new atmospheres and the yield 


point returns (strain ageing). 'Vhis type of hardening, caused by the assembly 


of mobile solute atoms round dislocations, is to be distinguished from the type 


of hardening considered previously in which dislocations are impeded by the 
strain fields of immobile irregularities, such as particles of precipitate. It is 
to be expected that, under suitable conditions, both types of hardening will be 
present simultaneously in the same material. 

The yield phenomenon is particularly marked in mild steel and soft iron, 
where it is known to be associated with the presence of small amounts of carbon 
or nitrogen (Edwards, Phillips and Jones 1940, Snoek 1941b, Low and Gensamer 
1944). Accordingly, it has been suggested (Cottrell 1948, Nabarro 1948) that 
the yield point of iron is due to carbon or nitrogen atmospheres. In the present 
paper this theory is examined by attempting to estimate the yield point and rate 
of strain ageing associated with carbon atmospheres; carbon is considered rather 
than nitrogen because it is more commonly present in iron, but the following 
discussion can be applied with very little modification to the case of nitrogen. 


§2. CARBON ATMOSPHERES IN «-IRON 


X-ray studies of martensite show that when carbon dissolves interstitially 
in #-iron it causes a marked lattice expansion and distorts the structure from cubic 
to tetragonal. ‘These large strains must lead to a strong interaction between a 
carbon atom and a dislocation, which is consistent with the pronounced yield 
point of iron. In the case of a (positive) edge dislocation a carbon atom can 
relieve hydrostatic stresses by entering the expanded region below the dislocation 
centre; also, as was first pointed out by Nabarro (1948), an additional interaction 
can be expected because the non-symmetrical character of the distortion produced 
by carbon should enable shear stresses to be relieved. For a screw dislocation, 
Frank (1948) has suggested that the interaction must be due solely to the relief 
of shear stresses because hydrostatic stresses do not exist round a pure screw 
dislocation, so that in this case atmosphere formation can only occur with those 
solutes which, like carbon, produce a non-symmetrical distortion. The conditions 
for a given solute to be able to produce the yield phenomenon may thus be more 
stringent than was previously suggested, and a high mobility and a large degree 
of misfit, although necessary, may not in themselves be sufficient. However, 
considering the atomic structure in a screw dislocation in an actual metal lattice, 
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_ it seems likely that the stress field will contain a small hydrostatic component, 
(so that a weak interaction may exist with atoms which cause only spherically 
iy symmetrical strains. 

a | The interstitial nature of the solution of carbon in iron and the non-symmetrical 
| distortion prevent the use of equation (1) for calculating the interaction energy. 


1) Since the axis of the tetragonal distortion produced by a-carbon atom is confined 


'© to one of the < 100 > directions of the body-centred cubic cell and a dislocation 
'# is confined to a slip-plane, with its direction of displacement lying in a slip 
i) direction, the system is not free to arrange that the tetragonal axis coincides with 


4) the direction of greatest tension in the dislocation field.. For a complete analysis, 


_ therefore, the orientation and form of a dislocation in the iron lattice must be deter- 
mined. The problem is much simplified if only the hydrostatic interaction is 
_ considered because hydrostatic stress is a scalar and the question of orientation 
4 is then avoided. In order merely to examine whether the anchoring force is 
“ sufficient to account for the yield point this simplification is justified and, accord- 


® ingly, the present treatment deals only with the hydrostatic interaction with an 
Sly P y Me 


edge dislocation. It is hoped to discuss the shear interaction in a later paper. 
a The interaction energy is determined by evaluating the work done by forces 
) from the dislocation field when a lattice expansion occurs equivalent to the 
) expansion caused by the introduction of a solute atom (Cottrell 1948). Assuming 
4 that each carbon atom alters the edges of the cell in which it is from a, a, a, to 
+ ©, a, a,, the mean changes in the lattice parameter produced by dissolving n 
8 carbon atoms in a crystal of N iron atoms are given by Ac=2(c—a)a/N and 
Aa=2(a,—a)n/N. Using the lattice parameter measurements on martensite 
) of Lipson and Parker (1944) and extrapolating to n/N =0-5, the volume change 
caused by the introduction of a carbon atom is given as Av=0-78 x 10-8 cm’. 
Hence, if the stresses due to the dislocation in the surrounding medium remain 
constant when this local volume change occurs, the interaction energy is 


ViZPNG == GO AUN an lumi (sae ee (3) 


| where p and © are the local pressure and dilatation, respectively, caused by the 
dislocation, and g is the bulk modulus. Using Koehler’s (1941) formulae for 
the stresses round a positive edge dislocation 


A 1—2v sine 
Ne lg re? 
and assuming that © =0 when r<7y we have 
Pe ay dy sins To. Vi=0 =F. 
r 37 l—-v 4 


The assumption that the dislocation stresses, and hence ©, remain unchanged 
when the carbon atom is introduced implies that the stress field is effectively 
constant over regions large compared with the size of the carbon atom, a reasonable 
approximation in regions remote from the dislocation centre. Near the centre, 
however, the dislocation stresses are rapidly varying functions of position and 
must be greatly altered by the introduction of the carbon atom. Tn the part 
of the crystal below the dislocation (i.e. where © is positive), which is the region 
of interest for the case of carbon in iron, the (negative) pressure must relax as 


the local volume change Av takes place so that the interaction energy will be less 
Vs 
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than is given by equation (4). Very near the dislocation centre equation (405 
will grossly overestimate the interaction and can only be used to give a rough | 
indication of the general order of magnitude. 

Taking G=8-28 x10" dyne.cm-?, A=2r,=2:-48x10-§ cm., and v=0-28) 
the constant A in equation (4) is 3x 10-?°dyne.cm?. At room temperature 
A/kT =76a. and, with such a value, the Maxwell—Boltzmann formula (2) gives 
impossibly large values of 2/m,) in the lower half crystal within about 10a. of the 
dislocation centre. The concentration in this region must be limited by other 
factors. An assembly of carbon atoms can never completely relieve all the 
dislocation stresses, but, as the atmosphere builds up, the progressive relief 
of stress will lead to a condition of saturation in which the addition of more atoms 
to the atmosphere no longer reduces the energy. An estimate of the number 
of atoms required for saturation can be made by equating the volume changes — 
produced by the atoms and by the dislocation field. ‘Taking a section consisting 
of one atomic plane threaded by the dislocation line, the radius R of a semi- 
circular annulus, in which the total volume change produced by the dislocation 
is equal to that produced by introducing one carbon atom, is determined by 


R 27 
Av=A| | Grete : 


Using the above value for Av and assuming 77 ~2a., Ris about 10a. It is thus 
reasonable to assume that in the equilibrium state only one carbon atom per — 
atom plane is contained within 10a. of the dislocation centre. Since the greatest 
tensile strain must occur immediately below the dislocation centre and in the 
atomic plane which constitutes the lower side of the slip plane, the most favourable 
position for the central carbon atom will be at « =37/2 and r~2a., and it seems 
likely that the central part of the atmosphere consists of a line of carbon atoms, 
parallel to the dislocation line, in this position. In the annulus between 10-20 a. 
the total volume change is again about the same as that produced by one carbon 
atom per atom plane, but this volume change is spread more uniformly over a 
much larger region of the crystal and, in this case, the carbon atoms could hardly 
be arranged along a line. ‘The atmosphere is pictured, therefore, as a central 
row of carbon atoms, situated just below the dislocation centre, which is surrounded 
in the lower half-crystal by a dilute distribution of the Maxwell—Boltzmann type. 
It should be noticed that an extremely small amount of carbon is required 
to provide atmospheres for all dislocations in a crystal; the actual value depends, 


Table 1 
Deno of dislocation (lines/cm?) see One 108 
% of C by weight absorbed f Full miei point 1032 105° 
in dislocations for No yield point Bs al One << Wiss 


of course, on the density of dislocations. ‘Table 1 shows the results of the assump- 
tion that a full yield point is obtained when every dislocation has one carbon 
atom per atom plane and that the yield point is absent at below one-tenth of this 
amount of carbon. 


$3) FORCE REQUIRED TO PULLEA DISTOCARIONSE ROM 
ITS ATMOSPHERE 
It is clear that the force required to pull a dislocation from its atmosphere 
is determined almost completely by the interaction with the carbon atoms of 
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4 ) the central row, since (i) these atoms, being so near to the dislocation, interact 
@ strongly with it, and (ii) the restraining force due to their interaction must reach 
. i its maximum at a very small displacement of the dislocation from its equilibrium 
‘ position, and at such a small displacement the contribution to the restraining 
1@ force from carbon atoms further out will still be comparatively small. From the 
@ point of view of calculating the yield point this is unfortunate, since it means that 
@) the important contribution of force is from atoms near the dislocation centre, 
‘y where the interaction energy cannot be accurately evaluated, so that only a rough 
| estimate is possible. 

__ In Figure 1 consider a positive edge dislocation lying along the Z axis (the 
X axis is the slip direction) and a line of carbon atoms CC’, lying parallel to the 
dislocation, in the most favourable position at a distance p below it. The inter- 
} action energy per atom plane, when the dislocation has been displaced a distance « 
@ in the slip direction is 


V(x)=—Asind/r=+Apfx®+p2)° EL. (5) 


and the corresponding force is F(x)=0V/0x=2Apx/(x?+p?)2 which has a 
maximum of 34/3A/8p? dynes per atom plane at x=p/4/3. 

The critical force for releasing a dislocation is thus Fy = 34/3A/8Ap? dynes/cm. 
and, by using the relation t= F/A (Mott and Nabarro 1948) where 7 is resolved 
shear stress, this can be converted to an equivalent critical shear stress 
T) = 34/3A/8A"p2 dyne.cm™ on the slip plane. Assuming that in the case 
of a tensile test the active slip planes are inclined at approximately 7/4 to the 
tensile axis, the critical tensile stress is 


SON OA DE AYNe.CIVaee | es ene (6) 


Table 2 gives the observed values of the yield point of annealed ingot iron 
at various temperatures, as derived from the graphical results of McAdam and 
Mebs (1943). The yield shows a strong temperature dependence and it is clear 


Table 2 
Temperature (° K.) ss = 85 144 195 303 373 
Yield point o (dynes/cm? x 10%) .. 6:07 4-14 2:83 1-45 1-10 
iin arc = om a 0-485 0-330 0-225 0-115 0-090 


that op» in equation (6) must correspond with the zero point yield stress. If a 
curve is fitted to these results and extrapolated to 0°K., the experimental value 
of a, is about 1-25 x 101° dynes/cm?. Using the precious values of A and A 
the theoretical value of co) can be made to agree with this if p is taken to be 74., 
which shows that the restraining force is strong enough to account for the yield 
point. If p is taken to be 2a., as suggested in §2, the estimated value of Gp IS 
an order of magnitude higher; however, it must be remembered that the inter- 
action energy is seriously over-estimated near the dislocation centre and the 
theory can only give a rough indication of the yield stress. Comparison with 
the experimental value of oy suggests that a reasonable value of A near the disloca- 
tion centre is 3 x 10-71 dyne.cm?. * 

Figure 2 shows the variation of interaction energy and stress with position 
of the dislocation for the cases p=2a., A=3 x 10-®° dyne.cm’, and p=2a,., 
A= 3 x 10-7! dyne.cm?. 
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§4. EFFECT OF THERMAL FLUCTUATIONS ON THE YIELD POINTE 


At temperatures above 0°x., thermal fluctuations help the external forces 

to tear a dislocation from its atmosphere and yielding occurs before o reaches oy. 
When a stress o = po, is applied the dislocation moves forward to the position x, 
of stable equilibrium (Figure 2) and can only be released by this stress after it 
has been taken further forward to the position x, of unstable equilibrium, beyond 
which the external force exceeds the restraining force. : 
The activation energy per atom plane for the release, by thermal fluctuations, 

of a dislocation which remains straight is thus the difference V(x.) —V(x,) in 


(2) «) 


o/0, 


ac in Angstroms. 
(b) 


Vac) in eV per Atom Plane 


Figure 3. Figure 2. Variation of interaction energy (_V) and applied: 


stress (o/o9) with dislocation position. 
Scale (1) A=3 x 10-*° erg.cm.: p=2 a. 
Scale (2) A=3 x 10-*! erg.cm. p=2 4. 


interaction energy minus the work done by the external force when the dislocation 
moves from x, to xX». In calculating the activation energy which thermal fluctua- 
tions must supply in order to release the whole dislocation, however, an energy 
contribution caused by the change in form of the dislocation must be considered. 

Let the plane of Figure 3(a) represent the slip plane. When the applied 
stress is zero the equilibrium position of the dislocation is ZZ’, vertically above 
the line of the atmosphere. Under the stress o the dislocation moves forward 
to the line x,. To produce yielding it is not necessary for the fluctuation to 
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| move the entire dislocation to x, as this would require an unnecessarily large . 
jenergy. If a loop ADC of the dislocation is thrown forward, and this loop is 
sufficiently large, the force pulling the looped part forward can exceed the 
» restraining force and the dislocation will then be torn from its atmosphere, the 
9 centre of the loop moving forward and the sides moving laterally towards the 
y= ends of the dislocation. This process depends upon the ability of the dislocation 
© to be bent and extended. The theory of flexible dislocations has been developed 
by Mott and Nabarro (1948) who have shown that a dislocation behaves like a 
™ one-dimensional soap film; it possesses an energy WL,.where L is its length 
) and W=1-—Sev. per atom plane, which it tries to minimize by becoming as short 
4 as possible. Thus to form the loop ADC, an activation energy has to be supplied 
/ which is made up of (i) the energy of interaction with the atmosphere, (ii) the 
# length energy associated with the increase in length, ADC — ABC, and (iii) the 
© energy gain due to the work done by the external forces. 

To calculate the exact size and shape of the smallest stable loop is difficult, 
© but it is clear that an essential condition is that part of the dislocation has to be 
brought by the fluctuation to at least the line x,, since otherwise there will be no 
part of the looped dislocation on which the net force acts in the forward direction. 
) Thus consider the simpler model shown in Figure 3(b). Here a triangular 
loop is formed, of depth BD, reaching to the line x, at its apex, and the value of 
the angle % is chosen to make the activation energy a minimum. ‘This loop is, 
) of course, too small to be stable; however, the true activation energy cannot be 
__yery much greater than is given by this model, for calculation shows that about 
one-half of the energy of the triangular loop arises from the length energy term 
and much larger loops could be formed with only a small additional increase 
in the length of the dislocation. 

It is assumed that the interaction energy of an element at any point on the 
portions AD and DC is the same as that of an element with that point as centre 
and lying parallel to ZZ’; this assumption is consistent with the results of the 
analysis by Mott and Nabarro (1948) of the strains round a zigzag dislocation. 
Accordingly, when an element dz of the dislocation moves forward from x, 
to x, with the formation of the loop, its interaction energy changes by 


(V(x) —V(x,)} dz = — cot W{ V(x) — V(x,)} ax 


and the contribution of the interaction with the atmosphere to the activation 
energy is 


2] —coty{V(x)—V(a)}dx. - 


~ Bo 


Taking V(x) = — Ap/A(x? + p?) erg/cm. this becomes 


a l= {tan-1(«,/p) —tan-1(«,/p)} + 2(x2 — x) Vr) | COUM MNS eigenen: (7) 


The increase in length energy is 
W(ADC — ABC) =2W(x,—x,)(cosecys—coty) ss eee (8) 


and, if the external force per unit length on the dislocation is F, the work it does 
when the loop is formed is 


UN Re COL S| uty.) Anew sits (9) 
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- Subtracting (9) from the sum of (7) and (8) gives the activation energy, U, which 
is of the form « cosecys—f cots; minimizing U with respect to % defines ¢ as 
cos1(B/x) and substituting this value then gives the activation energy in terms 
of Wand a/o,. It can be readily shown that U is of the form 

U =D(o/0)[AE(o/09){2Wp — AE(a/05)}]’, 

where D and E do not depend on the physical parameters 4, W and p. ‘Thus, 
if 2Wp> AE, U/(2AWp)! is essentially a function of o/o, only, so that altering 
A, W and p only alters the curve of U against o/o, by a scale factor, the form of 
the curve remaining unchanged. For the range of values of A, W, p and a/o, 
considered, AF is never greater than about Wp/2 and in most cases is much less 
than this, so that the form of the activation energy curve is almost indepéndent 
of the parameters. Figure 4 shows the curve for various values of A, W and p. 
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Figure 4. Activation energy for the formation of a loop. 
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It is to be expected that, as the stress is increased towards ap, the first successful 
attempt to form a loop occurs on a dislocation in a highly stressed region of the 
specimen. ‘The first stable loop to be formed is subsequently extended by the 
applied force until the entire dislocation to which it belongs is released from its 
atmosphere. This dislocation accelerates and enables other dislocations to 
break away by setting up an elastic disturbance in their midst; in this way yielding 
develops catastrophically throughout the entire specimen. This description 
agrees with the well-known observations that yielding starts in highly stressed 
regions and is propagated, in the form of Luder’s bands, from these regions 
along thespecimen. The elementary yielding process—the release of a dislocation 
from its atmosphere—occurs in both the start and the propagation of yielding; 
the difference between the upper and lower yield points is that, in the case of the 
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lower yield, the release of dislocations (at the edges of a Luder’s band) is helped 
by the elastic disturbance caused by nearby plastic flow, whereas in the upper 
yield this hi |p is not available. 


With the above picture the effects of temperature and testing rate on the yield 


point can be examined. When the stress in the neighbourhood of a dislocation 


lies in the range 0/5 to (a +do)/o, the time t, before a successful fluctuation 

occurs 1s proportional to exp{U(c/o,)/kT}. On the other hand, the time ¢, 

© spent in this stress range is inversely proportional to the rate of stressing, do/dt. 
_ Yielding should occur when ts ~t,, i.e. when the quantity 


S=(do/dt)exp{—U(c/o))/RT} sc aes (10) 

_ reaches a characteristic, fixed value. Thus, in experiments at constant testing 
me tate, the yield point should vary with temperature in such a manner that U/kT 
_ remains constant, and similarly, at constant temperature U should depend 
_ logarithmically on the testing rate. From McAdam and Mebs’ results (Table 2) 
the temperature dependence can be examined. If the experimental and theoreti- 
cal values of c/o) are equated at one arbitrarily chosen temperature, the value 
of U and hence U/kT can be found from Figure 4. Keeping U/RkT constant 
the theoretical o/c) values can be obtained at various temperatures from Figure 4 
and compared with the experimental values. Figure 5 shows the results obtained 


10 


— Theoretical 
o Experimental (McAdam & Mebs) 


0:8 


0-2 


(0) 100 200 300 400 
HK 
Figure 5. Variation of yield point with temperature. 
The theoretical curve is fitted to the experimental value of o/o9 at 195° kK. 


when the o/c, values are equated at 195°k. and the activation energy curve 4 
of Figure 4 is used. Since the activation energy curves have almost identical 
forms essentially the same results are obtained when the other curves are used. 
The good fit obtained over a wide range shows that the activation energy curves 


Table 3 
Activation energy curve oe 1 2 3 4 
W (ev. per atom plane). . sits 5 5 il 1 
Al (dyne.cm”?x10—*°). .. af 3 0-3 3 0-3 
OHS ee ae ae so Send) 107 140 47 


have the correct shape to account for the experimental observations. ‘The values 
of U/kT given by equating the o/o, values are shown in Table.3: Since Uj/RT 
in practice could hardly exceed about 50, it appears that the most reasonable 
values of W and A are those of curve 4. 
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In the present theory, the large temperature dependence of the yield point 
of iron is explained by the narrowness of the potential energy trough associated 
with an atmosphere of solute atoms (cf. Figure 2), so that only a small displace- 
ment of the dislocation, and hence a small activation energy, is required to overcome 
the restraining force. In a precipitation-hardened alloy, on the other hand, 
the existence of wide potential energy troughs causes a large activation energy 
for moving a dislocation, thus giving a small temperature dependence of yield 
strength (Mott and Nabarro 1948). 


§5. RATE OF FORMATION OF THE ATMOSPHERE 


In the present theory, strain ageing is regarded as the process whereby carbon. 
atoms in solution migrate towards free dislocations, at rest in an overstrained 
crystal, to form new atmospheres. To examine this process consider a free 
dislocation at rest in a region where the solute is initially randomly dispersed, 
with concentration m) atoms per cm*. Moving by thermal agitation the solute 
atoms acquire a drift velocity relative to the dislocation (Cottrell 1948) given by 


ve (DIETIVY @ A SaS Sa e (11) 


where D is the diffusion coefficient and V is given by equation (4). As the 
atmosphere builds up to equilibrium, both saturation and concentration differences 
become important in modifying the flow velocity, but in the early stages of atmo-. 
sphere formation they are negligible. As it is extremely difficult to take account 
of these factors quantitatively, only the initial stages of ageing, where the drift 
flow-is dominant, will be considered. 

Equation (11) shows that the atoms move normally to the equipotential 
surfaces of V. It is thus convenient to choose as coordinate lines for the problem 
an orthogonal, curvilinear set consisting of the families made up from the lines. 
of constant V and the lines of flow; the problem is, of course, two-dimensional. 
Since V=A sina/r, the first coordinate may be defined by the parameter. 
n =1/sin a, so that the equipotential lines are those on which 7 is constant. The 
flow lines, conjugate to these, are defined by the coordinate €=r/cos«. Using 
the standard transformation formulae, the elements of length along the € and » 
lines are, respectively, 


RRP eke. 7: (12) 
so 2p uy) 24a? 0.6 ere we 
an. 1S 
e2 4 2 
—A Bae CUNY palin 28) Grape ities: (13) 


where i is the unit vector in the 7 direction. As shown in Figure 6, the &, 7 
coordinate system consists of circles centred on the X, Y axes and passing through 
the origin. 

A carbon atom will migrate towards the most favourable position below the. 
dislocation, moving along a line of constant € and eventually come to rest at a 
position where its 7 coordinate has a value, — 7, of afew Angstroms. From (11),. 
(12) and (13), the time required for an atom to move an elementary distance in 
the 7 direction is 
_ RT Eb? 


dt = AD (@ +?) +72 dn, 
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1 and hence the time to move from (£, +79) to (&, —19) is 
PS lore lee @ Dnl(GXs tye lume (8 2 Noe 
= ape (|, oranlet+net + |” ert dnl(e +n99}) 
5 =(AT§*/AD)[m/2—tan- (no/£) + 3 sin {2 tan“!(q9/€)}]. 
| If 7 /€<1 this can be approximated to 
| t=nkT®/2AD 


@ since tan~'(n9/€) -~3 sin {2 tan“'(y9/E)} Sy o/E. After the time ¢, all lines of flow 
4 for which &<2ADt/zkT will have been drained of carbon atoms and will take 
) no further part in the ageing process. ‘Taking the smallest value of compatible 

with the above approximation to be 10., and using the value of §2 for A/RT, 


0 50 150 200: 
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Figure 6. Figure 7. Strain ageing of a low-carbon steel. 
(From results of Davenport and Bain.) 
©O =Ageing at 20° c. 
@ =Ageing at 40° c. 
wus — F}quation (16) assuming a carbon concen- 
tration of 0-003 % by weight. 


it can be seen that this conclusion should be valid for ageing times such that 
Dt>2 x 10-% cm®. As Figure 7 shows, strain ageing does not become appreciable 
until much longer ageing times than this have elapsed. 

Consider flow along a filament of width dé which is still supplying atoms at 
the time ¢. The time for an atom to move an elementary distance along this. 
filament, in the region where 7 = — 7, is 

RD oe? 
an a Le 
AD (P+ ai | 
and the area of the corresponding element of the filament is 
cna ADdé ; 
= In dé = dt. 
Cine are 
Using equation (14) to define the limits of integration, the total area which can 
supply carbon atoms to the dislocation in the interval of time from ¢ to ¢+d¢ is. 


aS —(2A Dt[nkT)* t (/AD\ dt 
all e-2 de + | dé) ie G) (=r) aoc (15) 


+(2A Dt/akTy> 


dt 
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An important feature is that, in the region of the dislocation field where the 


relation V =A sina/r holds, and outside the regions defined by equation (14) — 


from which the carbon atoms have been removed, the concentration of carbon 


remains constant during the drift flow. The hydrodynamical equation of — 


continuity, V.(pv) = —@p/dt, here takes the form V. {n(D/RT)VV}=0n/ot where 
n=n(r,«;t) is the concentration of solute. Since »=m)=constant at t=0, 
and since V2V =(d2/0r? +d/rdr + 0?/r°0a?)(A sina/r)=0 then dn/ot=0 and the 
condition =p is perpetuated. 

If unit length of the dislocation is considered, the number of carbon atoms 
which arrive within a time ¢ is, from equation (15), 


RNP LAD i Bees a\* (ADt\? 
weo=n2 (5) Car) Je 'et=nes(5) (Ger) 


~The total number, N,, of carbon atoms per unit length of the dislocation required 
to form an atmosphere of one atom per atom plane is 1/A, so that the degree of 
formation of the atmosphere at time ¢ is 


N(t)|N,=n)3(7/2)(ADERT). nee (16) 


§6. RATE OF STRAIN AGEING 


An experimental study of strain ageing has been made by Davenport and 
Bain (1935) who measured the increase in hardness of a low-carbon steel with 
‘time of ageing at various temperatures. Nabarro (1948) has shown that several 
features of strain ageing, shown by their results, are explicable on the basis of 
‘the present theory and provide evidence in its support. ‘To compare the ageing 
rates found by Davenport and Bain with the rate predicted in the previous 
-section the quantity H =(H,—#H))/(H,,—H >) which is a measure of the degree 
-of strain age-hardening, is first derived from the graphical results given in their 
Figure 11 (upper block); in this expression, H, is the hardness after ageing for 
-a time ¢, Hy is the value before ageing and H,,, is the maximum produced by 
prolonged ageing. ‘l’o assemble results from different temperatures in the same 
diagram it is convenient to plot H against Dt rather than the time ¢ itself, where 
D=5:2 x 10-4 € 9909/7 (Snoek 1941a, Polder 1945) is the diffusion coefficient of 
-carbon in «-iron, since all the results then lie on the same curve. _ 

Figure 7 shows the experimental and theoretical variation of H with Dt. 
“The theoretical curve has been derived from equation (16) by assuming that 
H= N(t)/N,, by using the values of A and A/kT given in §2, and by taking a value 
of my) corresponding to 0-003°% carbon by weight. ‘The assumption H = N(t)/N, 
-cannot be strictly valid since the first atoms to arrive at a dislocation ought to 
be more effective in anchoring it than those that arrive later. Equation (16) 
represents the observed ageing behaviour fairly well up to about N(t)/N,=0-3 
but beyond this range it gives a rate of ageing which is too high; this is to be 
-expected since the neglect of the effects of saturation and concentration differences, 
and the assumption H=N(t)/N,, will each cause the theoretically estimated 
rate to be too high in the later stages of ageing. 

The carbon concentration determined above by fitting the observed and 
-calculated ageing rates ought to correspond to the value for the residual carbon 
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in solution in «-iron at room temperature. Most of the carbon in irons and steels. 
of higher carbon content is, of course, precipitated as cementite and only that 
which remains in solution takes part in strain ageing. Experimental evidence 
for the solubility of carbon in iron at room temperature is meagre, but suggests 
that the above value of 0-003% is reasonable. Thus Yensen (1924) and Dijkstra 


) (1947) have given the upper limit of the residual solubility as 0:006°% and 0-001 °% 
respectively. 


§7, OTHER AGEING EFFECTS 


_ The migration of solute atoms into dislocations, causing the latter to become- 
anchored, ought to produce several mechanical and physical effects, which can 
_ be divided into two groups: (i) those resulting from the anchoring of the disioca-- 
tions, and (i1) those resulting from the removal of solute atoms from random 
solution to form atmospheres. 

The most striking effect belonging to the first group is, of course, the return 
of the yield point. In addition, a reduction of certain “‘anelastic”’ effects (Zener 
1946) is to be expected, and this is confirmed by experiment. A freshly strained 
spécimen possesses a high internal friction and shows slight plasticity at stresses 
well below the stress at which extensive plastic flow setsin. Zener has attributed 
these anelastic effects to local readjustments in position of dislocations in the 
_ slip bands of the worked material; this is reasonable since small displacements. 
__ of free dislocations in a Taylor (1934) array can occur under small stresses, even 
_ though a high stress is required to overcome the interaction forces of neighbouring » 
_ dislocations and so produce extensive flow. Low-temperature annealing, for 
times which would allow strain ageing to occur, greatly reduces this source of 
internal friction and causes a true elastic range to return, which is consistent 
with the suggestion that dislocations become anchored by the ageing process. 

Ageing effects belonging to the second group should be provided by properties . 
which are sensitive to the amount of carbon in random solution in «-iron. ‘Thus 
we may expect both electrical resistance and magnetic ageing effects. Recent 
work by Mr. A. T. Churchman and one of us (A. H. C.) has shown that a small’ 
decrease in electrical resistance of the right order of magnitude does in fact occur 
during strain ageing. Another effect has been observed by Snoek (1941a), 
who found that cold working and low-temperature annealing caused the elimina- 
tion of a certain type of elastic after-effect known to be caused by mobile carbon 
(or nitrogen) atoms in solution in «-iron. Snoek suggested that this treatment 
probably caused the solute atoms to become bound by internal stress fields 
created by the cold working. If we interpret the sources of the stress fields as 
dislocations, this suggestion becomes quite specific, since then the solute atoms 
cannot produce internal friction by Snoek’s mechanism because they are assembled 
in dislocations. 
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ABSTRACT. ‘The binding energy of a heavy nucleus is one of the quantities which a 
satisfactory nuclear force theory ought to predict correctly. Approximate calculations of 
this binding energy are performed, taking as nuclear force the static interaction of Moller 
and Rosenfeld. ‘The parameters appearing in the interaction have been fixed by data on 
light nuclei. 

The Fermi-gas model is taken as basis for the calculations, a first solution being obtained 
by a first-order perturbation (or variation) method. For refinement, asecond-order perturb- 
ation calculation is made, and Svartholm’s variation-iteration method is attempted. 
The first-order results yield about 10°% of the required energy, and the second-order about 

-40°%; a reasonable prediction of the size of the nucleus is obtained. The reliability of the 
perturbation results is uncertain ; Svartholm’s method does not appear to be well suited to 
the investigation of heavy nuclei. 

Two alternative modifications are made to the nuclear force to make it yield the correct 
binding energy according to the perturbation calculation. 


Sol SeLN RO DUC hin 


compatibility with the observed binding energies of heavy nuclei. By a 

semi-empirical procedure these can be approximately expressed by a simple 
function of atomic number Z and mass number A, containing a few parameters. 
One of these parameters, the coefficient of an energy contribution simply pro- 
portional to A, can be taken as representing the major effect of the nuclear forces, 
and a theory might be tested by the value which it predicts for this parameter. Of 
the theories which have been put forward, that of Moller and Rosenfeld (Moller 
and Rosenfeld 1940, Rosenfeld 1945, 1948) is capable of accounting for many of 


iL principle, one of the most important tests of any proposed nuclear force is its 


€ 
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the properties of light nuclei, and is at the same time relatively easy to handle: 
| this is the theory here considered in relation to heavy nuclei. The parameters 
__ appearing in the interaction are fixed by data on light nuclei. The problem posed 
__ is the solution of the wave equation for the ground state of a very heavy nucleus. 
Approximate methods must be employed for this, and in consequence the re- 
_ liability of the test cannot very well be assessed. 

The methods used are based on the Fermi gas model (§ 3.1), a first approximate 
solution being obtained by a first-order perturbation (or variation) treatment 
) ($3.2). Two different possibilities of refinement have then been attempted, 

__viz. a second order perturbation calculation (§ 3.3), similar to that first performed 
_ by Euler (1937) for a different form of interaction; and Svartholm’s variation- 
| iteration method (§5). The former leads to positive results, but it appears that 
Svartholm’s method is not well suited to the investigation of heavy nuclei. The 
results of the perturbation calculations are presented in § 2, and the calculations are 
criticized in §4. 
§2. FORMULATION OF PROBLEM AND RESULT 
2.1. Formulation of Problem 


The binding energy — £ of a nucleus of atomic number Z and mass number A 
can be represented approximately by the semi-empirical formula (Rosenfeld 


1948, Chap. 2) 
N-Z\? EN ek. 
B= —¢,A| 1-c4/ A ) |eeate 3(2)5, OOS (1) 


where ¢,, Cg, c, are constants, N is the number of neutrons, A—Z, and 74 is the radius 
of the sphere whose volume is that occupied by one nucleon, viz: 


79 =(radius of nucleus)/A?. 


7) is observed to be approximately constant, and is taken to be exactly so in 
equation (1). 

The several terms in the formula have the following significance: The first 
(if cp is neglected) represents the bulk of the “nuclear forces”’ binding the nucleus. 
The factor containing c, is a correction for disparity in the numbers of neutrons 
and protons. ‘The second term is a “‘surface-tension”’ effect, and the last term 
represents the electrical energy of repulsion of the protons. 

We are chiefly concerned with the parameter c,, whose observed value, found 
by fitting equation (1) to the nuclear energy surface, is to be compared with that 
predicted theoretically. By extrapolation of equation (1), c; may be interpreted 
as the binding energy per nucleon of an ‘‘ideal’’ nucleus, 1.e. an assembly of an 
infinite number of neutrons and protons in equal proportions, if the electrical, 
repulsion between the protons is removed. 7) can also be estimated theoretically 
and compared with the observed value. 

For the theoretical analysis the static interaction of Moller and Rosenfeld 
(Moller and Rosenfeld 1940, Rosenfeld 1945, 1948) is employed. At this point, 
dimensionless variables, to be used throughout, are introduced, by setting up the 


following ‘‘ absolute” units :— 


Unit of mass—mass M of a nucleon (proton and nev- | 


tron being assumed equally heavy). 
Unit of length—range of nuclear force, 1/A; defined by ( +++--- (2) 
A=pc/h where ~=mass of meson. 
Unit of angular momentum—/. 
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These yield the derived “absolute” unit : 


Unit of energy—Mc?(/M)? =A?2h?2/M. na ea (3) 
Then in these units, any pair of nucleons i and j has the mutual potential 
wv5i= (7, ° t)[g? +f(6; > o;)| exp ( - Ei) /&ij se icuatences (4) 


where a; and 7; are the spin and isotopic operators respectively for the i-th nucleon, — 
&,; is the separation of the nucleons, and g and f are interaction constants. 

The values of the parameters 1, f and g are taken from Frohlich et al. (1947, 
equation (76)), where they were chosen so as to fit correctly the two lowest energy 
levels of the neutron-proton system and the ground level of H°: 


b= 220 i.e Soe 27 = 2501 rt eg ce (5) 
These make the unit length, 1/A, equal to 1-75 x 10-% cm.; and the unit of energy 
equal to 14-4mM.M.u., or 13-4Mev. 
The observed values assumed for the quantities c, and ry of equation (1) are. 
those given by Bethe and Bacher (1936, § 30), viz. c, = 14-9 M.m.u. (1-03 in absolute 
units), 7) =1-48 x 10-4 cm. (0-85 in absolute units). 


2.2. Results 


The “‘ideal”’ nucleus is subjected in § 3.2 to a variation treatment, on the basis 
of the Fermi-gas model. This yields a value of the binding energy per nucleon 
for each initially arbitrarily assumed size of the nucleus. Minimization of the 
energy with respect to the size yields estimates of both the actual binding energy 
and the actual size. 


In Figure 1 are plotted the expectation values of the kinetic, potential and total 
AL 


| Figure 1. Energies per nucleon calculated 
@ on Fermi-gas model, vs 79, representing 
size of nucleus. 

(a) TIA; (6) WIA; (©) (T+W)A; 
(d) E’/A ; (e-) (T+W+E/A. 

The kinetic energy T/A (curve (a)) and. 
potential energy W/A (curve (6)) almost 
cancel. The maximum binding energy 
—(T+ W)/A by the variation (first order 
perturbation) method is far below the 
observed value on the radius 7) which 
minimizes (7+ W) is in good agreement 
Observed bindi with the observed value. 

The second order perturbation energy 
E”/A_(curve (d)), when added to 
(T+ W)/A (curve (e)), raises the binding 
energy considerably, but it is still well 
below the observed value for all 7. 
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energies (curves a,b, c, respectively) per nucleon as functions of 79, the “‘ radius per 
nucleon”. At the minimum of the total energy it is found that the kinetic energy 
per nucleon = 7/A =0-97, the potential energy per nucleon = W/A = — 1-06, and 
therefore total energy per nucleon= —0-09. (All energies are expressed in the 
absolute units of equation (3)). ‘The binding energy is thus of a smaller order of 
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} magnitude than the value 1:03 required, as in the case of similar calculations 
7 previously carried out with different forms of interaction (Bethe and Bacher 1936, 
| p. 157, Rosenfeld 1948). The minimizing radius ry is equal to 0-90 (in the 
# absolute unit of equations (2)), which compares well with the observed value 0-85. 
Although the binding energy obtained is so small, it may be noted that a quite 
f moderate increase in the strength of the potential would increase it rapidly, 
}) because the potential and kinetic energies under the assumed conditions almost 
cancel. 
' The second order perturbation calculation, which is carried out in $3.3, 


i yields a term E” td be added to the former energy (T+ W). The energy is again 
| found for variable size of the nucleus; although perturbation theory does not 


| strictly prescribe what size will give the best prediction for the energy (T+ W+E’), 
_ the minimum of the latter quantity with respect to size was taken as being probably 
+ asgoodasany. However, there seems no reason in principle why the correspond- 
| ing minimizing radius 7, yielded by the second order calculation should be any 
) better as an approximation to the actual radius than that yielded by the variation 
calculation. ‘The latter has therefore been retained as the theoretical prediction 
) of the radius 7. 
Graphs of E”/A and the total energy (+ W-+£")/A are shown in Figure 1, 
curves (d)and(e). It will be seen that| £”| is much larger than| 7+ W| at all radii, 
| indicating that the variation treatment is extremely poor. The minimal value of 
the total energy per nucleon is found to be—0:41. ‘This is much better than before, 
| but still falls far short of the required value —1-03. This second result, then, 
| indicates that the Moller-Rosenfeld interaction is too weak for the binding of 
| heavy nuclei, but it is not clear how much in error the calculation still may be. 
If the second order calculation is tentatively accepted as substantially correct, 
one may investigate what slight modifications to the nuclear force could yield the 
observed binding. ‘Two possibilities are considered. Firstly, the interaction 
— effective in heavy nuclei may differ from that in light, owing to ‘‘many-body”’ 
— forces. ‘To represent this effect the magnitude of the potential has been increased, 
the parameters f and g of equation (4) being varied in the same proportion, while 
_ the meson mass, and hence 1/A, are left unaltered. ‘The resulting energy per 
— nucleon in an “‘ideal”’ nucleus, taken to be the minimum of (+ W+£")/A with 
respect to variation of size of the nucleus, is plotted against f? in Figure 2, curve (a). 
It will be seen that f? requires to be raised from the assumed value of 0-595 to 
 0-73—approximately 25°/, increase—to yield the observed binding. The 
corresponding predicted size of the nucleus, given by the 7) which minimizes 
(T+ W), is plotted against f? in Figure 2, curve (6). When f* is raised sufficiently 
to yield the correct binding, 7, falls to 0-65, only about 75°% of the observed value 
0-85. Thus if many-body forces are invoked to account for the binding energy, 
_ the agreement regarding the size of the nucleus deteriorates. 

Secondly, the interaction can be altered so as to yield the correct binding energy 
per nucleon in an “ideal” nucleus, by making this one of the quantities by which 
the parameters of the interaction (equation (4)) are fixed. The interaction can 
then be compatible with only two at most of the three quantities pertaining to light 

~ nuclei, by which the parameters were fixed by Frohlich e¢ al. (1947): the two 

lowest levels of the neutron—proton system have been retained, while the ground- 

state energy of H? was abandoned. By fitting the two former quantities, f* and g* 
PROG, PIIYS. SOC. LXII, I—A S) 
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were calculated approximately for each value of the meson mass p, a method 


differing only slightly from that of Frohlich et al. (1947, §3) being employed. 
Then for each meson mass the binding energy per nucleon of an “ideal” nucleus 
was calculated, as shown in Figure 3, curve (a). Again, the minimum of 
(1 +W-+E") with respect to size of the nucleus was taken. It is found that to 
obtain the observed binding energy the mass of the meson required is 400m. 
The values of f? and g? appropriate to this meson mass are f? =0-60, g? =0-11. 
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Figure 2. Effect of postulated “‘ many- 


SR Bae ae Figure 3. Predictions for heavy nucleus, when 


meson mass is variable, and interaction para- 


energy and size of nucleus with 
effective nuclear-force parameter f?. 
Curve (a): Binding energy per nucleon 
by second order perturbation calcu- 
lation, i.e. (—1) X minimum of 


meters are chosen to fit two observed constants 
of neutron—proton system. 

Curve (a): Binding energy per nucleon by second 
order perturbation calculation, i.e. (—1) X mini- 


mum of (T+W+E’")/A. 


(T+W+E")/A. 
Curve (6): Radius 7) which minimizes 
the variation energy (T+ W). 

To obtain observed binding 
energy, f” has to be raised from 
given value 0-595 to 0-73: the 
radius 7 then falls to 0:65 (observed 
value 0°85) 


Curve (6): Radius 7) which minimizes the variation 


energy (T+ W). 

To obtain observed binding energy, mass of 
meson must be 400m: rp then falls to 67% of 
observed value. 


Together with each binding energy is found a prediction of the size of the 
nucleus (the value of 7) minimizing (T+ W)), which is shown plotted in Figure 3; 
curve (b). When the meson mass is 400m, 79 is only 67°, of the observed vel 
Thus, again, fitting of the binding energy is incompatible with correct nuclear 
size. 


CALCULATION OF THE, BINDING-ENERGY 


3.1. Use of the Fermi-gas Model 


On this model the nucleus is treated in the lowest approximation as an assembly 
of non-interacting particles, enclosed in a box at zero potential, with high potential 
walls. If the nucleons—neutrons and protons alike—are described by space, spin 
and isotopic coordinates, then the exclusion principle requires the wave function 
to be antisymmetric with respect to interchange of all the coordinates of any pair 
of particles. Suitable antisymmetric wave functions for the system may be built 
up in determinantal form from different individual wave functions for single 


§ 3. 
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\) particles. The composite wave functions so obtained will be approximate 
¥) eigenfunctions for the system, provided the individual wave functions are given by 


Y= jo, are eterele (6) 
#, 1s a function of the space coordinates only, zero outside the box, and 
ieeed eS CxD RSE) oe pe asec (7) 


inside the box. Here 2L is the length of the sides of the cube constituting the box, 


- € is the space-coordinate of the particle concerned, k, is the dimensionless mo- 


‘mentum, quantized by the condition 
: k,=n,7/L, safe (9) 


where n, is any vector whose components are all integers. w, in equation (6) is 
one of the four spin and isotopic wave functions: 


Oe = 26, Opel ee =O nO, =, eae (9) 
where « and f are the Pauli spin wave functions, and é and 7 are similar for the 


_ isotopic coordinates. 


To represent the ground state of the system, the composite wave function is 
built up of the A individual wave functions V’, of lowest momentum k,. For 
the ‘“‘ideal’’ nucleus considered, the A/4 spatial wave functions 7%, of lowest 
momentum are associated each in turn with all four w,,. 

The Fermi-gas model described has now to be used to calculate the ground- 


' state energy of the system of nucleons with nuclear interactions operative. 


In §3.2 a variation calculation is performed, using the ground-state wave 
function of the model, together with the Hamiltonian of the “ideal” nucleus, to 
find 

ems Py tee 1 tee) Ve (10) 


where T is the expectation kinetic energy, and W the expectation potential energy. 
This is carried out initially for arbitrary dimension L of the box. 

A perturbation approach is adopted in §3.3. Here the Fermi-gas of non- 
interacting particles is taken as an unperturbed system and the nuclear interaction 
W is applied as perturbation. The unperturbed energy is the kinetic energy T, 
and the first order perturbation energy is W, both already calculated. The 
second order approximation adds an energy EL”, given by 

Pe NW yn ta ols i eee (11) 
2 ; 


where the suffix R indicates the Rth unperturbed state and the suffix 0 the ground 


state. ’ indicates summation over all states but the ground state, obtained by 
building the composite wave function from all possible combinations of different 
individual wave functions. The size chosen for the box is again initially arbitrary. 
The index 0 on the energy indicates the unperturbed kinetic energy of the state 
concerned. 

3.2. Variation Method 


The expectation kinetic energy T is first calculated. In the dimensionless 


variables employed, it is found that 
A 
T= -450*( ys v5) Oe pe ae (12) 

=1 


s states 
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where ® is the ground-state wave function of the Fermi-gas model and fdq | 
indicates integration or summation over all coordinates. The individual spatial — 
states occupied form a lattice ink space, with a spacing 7/L (cf. equation (8)). The » 
occupied points, each four times occupied, are all those points falling within a 
sphere of radius a, given by 


_ At/3_,\t (9a) a 
a= 37(57) = 7 | a ee © ey Bokeromete te ( ) 


The sum in (12) can approximately be replaced by an integral in k space, evaluation 
of which yields 


T/A =3a?/10. ar ee (14) 


The expectation potential energy is next found, using for the potential operator 
= w;, where w, is given by (4), and the same wave function ®. ‘The result is 


Ei 


W=42 [(rs|w|7s)—(rs|w|sr)], == aeeeee (15) 
where, for instance, bei 
(rs| w| rs) = J PECL) VF (2)e012'¥',1) ¥5(2) dgn age 


and the sum takes r and s over all occupied individual-nucleon states. 


The matrix elements involved may be split into spatial and non-spatial parts” 
if we write (cf. (4)) 


Dj =(%;.7;)[2? +/7(4- 9;)] ; Uj = {exp (— Sic esecee (16) 
so that (15) becomes ; 
W=3 i{ z (ur| | Ww») | E (Im| o| im) | se | z (ur| P| vu) || Sm v| ml) |) 
: | _ se. (17) 


, v are non-spatial, and /, m are spatial quantum numbers (cf. equation (6)), the 
latter to be summed over the occupied sphere inkspace. The first, non-exchange, 
term vanishes, on account of the non-spatial factor, in conformity with the satur- 
ation requirements of nuclear forces (Rosenfeld 1948, chap. XI). ‘Towards the 
evaluation of the remaining exchange term, it is found that: 


4 - 
2 Gav] Tye) 12 (22263) eee (18) 
il 
and 
7 1 
 |@| wal) = ————-—-——. neers 19 
(lm \v| ml) FEAR SOT SG (19) 
The sum over /, m in (17) is replaced by a double integral in k space, leading to the 
final result 
W/A = —(9/87)(g2 + 3f2)Q,(a), 
where 


Q,(a) =2 E + soar | In (1 +(2a)2) + (2 - a) =F ime Qa) eee (20) 


The kinetic and potential energies of equations (14) and (20) are plotted against 
ry in Figure 1, curves (a) and (6), for the values of f? and g? given in equation (5). 
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3.3. Second Order Perturbation 


The second order contribution to the energy (equation (11)) requires evaluation 
of the matrix element W , for each excited, unperturbed state R. The only states 
giving non-zero contributions are those built up from a set of individual wave 
_ functions ‘’, and differing only by one or two constituents from the set constituting 

_ the ground state; and of these, the only ones which may not be neglected are those 
1, obtained from the ground state by replacing two individual wave functions ©’ Rot be, 
i, by two others ‘l’,., ‘V’, satisfying a “conservation of momentum” condition 


BB Ek keko + 2 ee wee (21) 


I] 
t 
i 
Be 
/ 
| 


_ where /, m etc. are the spatial quantum numbers pertaining to the individual states 


| 7, s, etc. respectively (cf. equation (6)). . Development of (11) then yields 
RO PX iio ’s’)*(rs|w|s’r’) 
ny yr Ls let's’) | (rs|w|r’s’)*(rs| Lay: 
y 1,87 58" (Ri + Riv) a (k; 1" Rin) 
Here & means double summation over all ground-state constituent individual 
states ¥, and ‘’,, and &’ means double summation over all the other individual 


TERE 
states ‘’,, and ¥’, conforming to the selection rule (21). The spatial and non- 
_ Spatial factors in the matrix elements are now separated out, as in the previous sub- 
section. The sum over the non-spatial quantum numbers involved is readily 
_ carried out, with the result 


ba =k , A\(dn || U'm’) 2 — B(Im| 0 | I'm’)*(lin| 0 |m'l’) - 
ee sso CeDaC rc arene | a 


where A=48(e4+3f4), B=12(3f4—6g?f?—g*). The spatial quantum numbers 
/, m are each to be summed over the sphere of radius a in k space; and /’, m’ are 
each to be summed over the space exterior to this sphere, subject to the selection 
rule (21). 

Evaluation of (dm|v|l'm’) yields 
7 1 
L3 1 +4|k,—k,—-k,,+k,,? \ 


(lm| v|l'm') = 3 


The sums in (23) may, as usual, be replaced by integrals in k space. By an 
obvious abbreviation, (23) may be written 


Pies AL OBI times ae, terry tgs (25) 


It has been found possible to evaluate J but not J. However, it appears that the 
intractable term containing J is negligible; for it readily follows from the form of 
the integrals that J>J>0, and insertion in (23) of the values of f? and g” given 
in (5) shows that | A|>|B|, in fact A=217, B=3-6. aie 

Thus neglect of the term BJ introduces an error of less than 2%. KC 
smallness of the error is to some extent due to fortunate given values of f? and g?: 
it may be noted, however, that the form of A and B is such that the error could not 
exceed —25°%, to +50°% for any values of f? and g* whatever. 

Evaluation of J from (23) and (24) then yields finally 


E"/A =(—9/2n?)(g4 + 37!) 224); 
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13. 4 1 7 le 5 
Qy=In(1 + (24)).| 75 - zin2+ aan |= 73 * 7g in2t¢ In2. Bayi 


oe I 11 1 7 semee 
aa. pe eves d 
+cot (2a)| 2a 5 oa |+ | 5 (2a) + 5 +5 +5a0 I, In (2. cos x) dx 


i 1Jn (1 —x?) af). sau x 
+2[ 145055 |[ aya 5 (2a)® | at aaa en (26) 


This is plotted against 7) in Figure 1, curve (d), for the f? and g? of equation (5). 


§4, CRITICISM OF (CALCULATIONS 


The argument has proceeded by the three following stages, each of which 
requires to be examined: (i) by a semi-empirical approach, based on equation (1), 
the binding energy per nucleon of an ‘‘ideal”’ nucleus has been determined; 
(ii) the parameters in the nuclear force employed have been chosen to fit certain 
properties of light nuclei; (iii) the nuclear force has been used to calculate approxt- 
imately the binding energy per nucleon of an “ideal” nucleus, for comparison with 
the result of (i). 

The procedure of the first stage is well known, and an extensive literature 
(e.g. Bethe and Bacher 1936, §30, Mattauch and Fliigge 1942, p. 90, Feenberg 
1947) establishes that the results obtained for the ‘ideal’ nucleus can be relied — 
upon to the accuracy relevant here. 

Stage (ii) was carried out by Froéhlich e¢ al. (1947), and is again probably 
sufficiently accurate as a basis for the present calculations. | 

Stage (iii) concerns the accuracy of the approximate solutions of the Schrédin- 
ger equation. ‘The inadequacy of the variation calculation is palpable, in view 
of the relative largeness of the second order perturbation result. As regards the 
latter, it may be noted that, even if high order perturbation energies were calcu- 
lated, and if they converged, the Fermi-gas model could only yield the energy of 
the system of particles when enclosed in a box of the assumed size, with infinitely 
high walls. ‘I"he method is thus best considered as an approximation to a solution 
for the nucleus so immured. For each order of perturbation, the best energy 
would be taken as the minimum with respect to variation of the size of the box. 
From this viewpoint, three orders of perturbation energy have been obtained, viz. 
T, Wand E”. ‘The only available criterion of the goodness of the approximation 
is the rate of apparent convergence of these terms; per nucleon, they are in order 
(for the interaction parameters assumed, and at the radius which minimizes their 
sum) 1:35, —1:41, —0-35 (sum = —0-41) (in absolute units). These figures might 
suggest inception of convergence, but only at a slow rate. 


§5. SVARTHOLM VARIATION-ITERATION METHOD 


The method developed by Svartholm (1945) has been used by him to investi- 
gate the ground states of the lightest nuclei, and its application to heavier nuclei 
seemed worth considering. It is found, however, that the method loses its usual ~ 
advantage of assured convergence in such an application; and in practice a 
calculation with the Fermi-gas model has yielded no useful results. 

The convergence question will be dealt with first. Reference should be made 
to Svartholm’s thesis (1945) for an account of the method, but it may be remarked 
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ng here that by an iteration procedure one may obtain in principle a succession of 


a4 approximations Noy A, .... to the eigenvalue A of the interaction constant required 
g to yield a given binding energy. A virtue of the method is that, subject to certain 
) conditions, the series of approximations Aj, A,.... decreases monotonically 
| towards the required eigenvalue. One of those conditions, however, is that the 
) potential function wu should be “‘negative definite”’, i.e. for any arbitrary wave 
| function ®, {®*u®dq<0. Now one requirement of a nuclear force suitable for 
# heavier nuclei is that it shall lead to the “saturation” phenomenon (Rosenfeld 
| 1948, chap. XI), and it can be shown very generally that this implies that the 
m potential shall not be negative definite. Thus for such a force, there is in general 
¥ no assurance of the convergence of the series of approximations. 

Nevertheless, the author has made a calculation with the method, using as 
initial wave function the Fermi-gas function of §3.1. The zero-order calculation 
yielding Ao, is then simply equivalent to the variation treatment of §3.2. A, was 
next calculated, only leading terms with respect to powers of A being retained. 
The result found was that, quite generally, for any type of saturating short-range 
potential, A, is identically the same as Ay. ‘Thus to this order, the Svartholm 


treatment gives no improvement over the variation treatment. 
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REVIEWS OF BOOKS 


Practical Five-figure Mathematical Tables, by C. Attwoov. Pp. v+74. 
(London: Macmillan and Co. Ltd., 1948.) 3s. 


_ The publication of this set of tables is an event of more importance than would at first 
meet the eye. The first mathematical tables were naturally of many figures, and then, in 
the nineteenth century, logarithmic and trigonometric tables settled down to seven figures 
as the norm. Towards the end of the century, and for the first part of the twentieth, the 
value of four-figure tables for the scientist was realized, and many a set was issued. With 
them, were a few six-figure sets and a fair number of five-figure ones. But their compilers 
followed each other rather monotonously, and seldom introduced real innovations, except 
for a few impracticable ones. The present reviewer, for example, has (or has lent and lost), 
more than a dozen different books of four or five figure compilations, most of them 
indistinguishable from each other. 

The difference that marks this, and that ought to make it very widely known and used, 
is that the author has thought out afresh what functions are wanted, at what intervals they 
should be tabulated, and the most convenient means for interpolation in them, 
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His choice of functions does not in fact offer any surprises. There are no transcendental 
functions such as we find in Dale’s five-figure tables, or in Jahnke and Emde, but a generous 
selection of those likely to be wanted in an all-purpose compilation. Logarithms and 
anti-logarithms, the six trigonometric functions, both natural and logarithmic, for degrees 
and minutes, we may take for granted. In addition, Mr. Attwood gives us logarithms of 
reciprocals (to save subtraction) and trigonometric functions of angles given in radians ; 
he also tabulates squares and square roots, cubes and cube roots, the fourth and fifth roots 
of integers to 100, and the exact values of powers up to the sixth of integers from 1 to 100, 
as well as reciprocals. ‘hen there are natural logarithms and exponential functions and 
(a concession to engineers ?) the areas of circles. 

Logarithms and trigonometrical functions offer difficulties in places, as we all know, 
where they change rapidly, and there are many ways in use for turning the difficulty. We may 
use the reciprocal or the function of the complement of the angle, or we may tabulate some 
function which varies less rapidly, we may vary the interval of the argument, or we may use 
the inverse function. Mr. Attwood points out, at the appropriate page, when the former 
expedients are desirable, but in general he relies on the method of changing the interval. 
This he has done most judiciously. Related to this matter is the question of the use of 
proportional parts : some five-figure tables give proportional parts in little marginal tables, 
leaving the user to form his own differences and to take the appropriate proportional parts 
(the reviewer’s favourite five-figure tables, due to Schlémilch, take them to one extra decimal 
place and indicate whether the tabular entries are in excess or defect, so that there is no loss 
of accuracy) whilst others give mean proportional parts applicable to a whole line of the table. 
These are dangerous when the differences change rapidly, and it is not uncommon to give 
two lines of mean proportional parts, splitting the tabular entries to share them between the 
two. Inthe present volume, this device is widely adopted and somewhat extended, so that 
as many as four lines of mean proportional parts (m.p.p.) may be given to a line of the table. 
The m.p.p. have been individually examined to get the best compromise, and are printed in 
red whenever they have to be subtracted. ‘lhe author asserts that more than 95°% of the 
160,533 combinations obtainable with the m.p.p. are in error by 0 or 1 unit, only 0:6% 
leading to errors of more than 3 units. Where there is a danger of error of five units, the 
m.p.p. are in italics, so that a user can use proper interpolation methods if needed. There is 
a table showing the maximum possible error in every table, but here the author is unfair to 
himself. No one could use m.p.p. to stretch across more than half the tabular interval, 
but should always take differences from the nearest entry. [his would reduce the errors in a 
number of places. Where m.p.p. are not given, the user is instructed to adopt linear or 
second-difference interpolation, and is given instruction in this, with a critical table of Bessel 
coefficients, but second differences are never printed as part of the tables. 

To close with a few words on general points. The author has varied in his interpretation 
of five figures. In some tables, it means five decimal places, as in the cosine of 89° 42’, 
which is 0:00524, and sometimes five significant figures, as in the secant of the same angle, 
which is 190-99. From the latter, with the aid of the table of reciprocals, we may find that 
the cosine is really 0-0052359, or from the tables of logarithmic cosines and antilogs, we may 
find it as 0:0052360. In the table of secants, there is a change at 45° from five decimal to 
five significant figures. 

The printing is good, and there is a strong limp cloth cover. If any change in a good 
book might be suggested, it is a thumb index to enable two tables to be used in conjunction 


more easily. To find out if there are misprints will call for longer use, or for deliberate and ~ 


systematic examination of differences. J. Howe 


Some Aspects of the Luminescence of Solids, by F. A. KrOGrEr. Pp. xii+310. 


(New York and Amsterdam: Elsevier Publishing Co., 1948.) $5.50 (30s.). 


The luminescence of solids is a subject which has benefited much from the development 
of theoretical conceptions of the solid state during the last fifteen years ; perhaps before 
long we shall see new theoretical ideas arising as a result of new experimental work. The 
work of the Lenard School during the early years of the century was the most significant 
in this field, but, unfortunately, the lack of a satisfactory theoretical picture of the behaviour 
of electrons in solids robbed much of the work of the value it might otherwise have had, 


ce gs, e+ a 
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During the last ten years a great deal of attention has been given to the study and under- 
standing of afterglow or phosphorescence phenomena in solids, and there have been notable 
_ advances. While to some extent a reasonably satisfactory picture can now be given of 
| these effects, we know little as yet about the nature of luminescence centres and electron 
traps in solids. These problems_in relation to infra-red sensitive phosphors and their 
| properties are particularly important, and it is in this direction that we may hope for new 
of advances. 

Dr. Kréger’s volume is of a rather special category, in that it is a research monograph 
and not a comprehensive review of a wide field of work. The volume is based on experi- 
mental work carried out in the physics laboratories of the N. V. Philips Gloeilampenfab- 
rieken, Eindhoven, during the last five years, or it may be assumed that it represents a 

-convenient form of publication of work carried out during, and held up by, the war. 
The first chapter contains a useful account of energy levels in pure and perturbed 
crystals, and an attempt is made to develop a general picture capable of covering all possible 
luminescent effects. Some of us will think that it is a little early for such attempts, but are 
} nevertheless grateful for the setting out of Dr. Kréger’s ideas in coherent form. 
7 In the following four chapters the author discusses new experimental results on 
_ tungstates, molybdates and phosphors activated by manganese, titanium and uranium. 
In the final chapter the effect of temperature on the efficiency of luminescence has been 
considered. ‘The data presented in these chapters will be of value to other workers in the 
subject, but can hardly be expected to interest a wider audience. 'The book ends with 
36 pages of tabular matter of the kind only too familiar to readers of Vol. XXIII/1 of 
the Handbuch der Experimental Physik. Is it of much scientific value to be told that 
- the luminescence of magnesium orthosilicate is ‘‘ violet’? ? Or that the luminescence 
is “strong” or “weak” or “‘faint’’? I think not: it is time for workers and writers 
in this field to rid themselves of these older inexact traditions if the subject is to advance 
as it should. Dr. Kréger’s volume is a useful addition to the research literature of the 
luminescence of solids, and the work described must, for the most part, have been carried. 
out under appallingly difficult conditions. J. T. RANDALL. 


- Electrons and Nuclear Counters, by 5S. A. Korrr. Pp. xi+212. (New York 
and Toronto: D. Van Nostrand Co. Inc., 1947; British Agents : 
Macmillan and Co. Ltd., London.) 18s. net. 


Until a few years ago, an aura of mystery surrounded the mode of operation of a Geiger- 
Miller counter and it was regarded as an erratic device which had to be made according 
to some favourite ‘‘ cookery-book”’ recipe. Recent advances in our knowledge of the 
discharge mechanism in these counters has done much to clarify the situation so that they 
are now standard tools for work on nuclear physics and are being satisfactorily manufactured 
commercially. There has been a deeply felt need for a book giving an up-to-date account of 
the state of the subject and Professor Korff is to be congratulated upon a pioneer effort on 
these lines. 

It should perhaps be stated at the outset that he has produced a survey suitable for the 
student and general reader rather than a work of reference for the specialist. The book 

~ opens with an introductory chapter describing the different modes of operation of an ion 
counter, that is, its use as an ionization chamber, proportional counter or Geiger-Miller 
counter. The next three chapters, which form the main body of the book, discuss each of 
these modes of operation in more detail, the main stress being laid on the Geiger-Miller 
region. The fundamental principles involved are well brought out and a good general 
account is given of the discharge which occurs in the gas of the counter and of the various 
factors which influence it. Unfortunately this account suffers from a rather poor arrange- 
ment of the material so that there is a considerable amount of repetition and frequent and 
very irritating references forward to matters dealt with subsequently. 

In chapter five, the author gives a very clear account of the principles and procedure 
involved in the practical construction of counters. He has dealt only with the cylindrical 
cathode, axial-wire anode type of counter. Sections dealing with other types, such as point 
counters and parallel plate ionization chambers would have been valuable additions. ‘There 
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follows a brief chapter on the errors and corrections involved in counting. The treatment 
is rather sketchy and can only be regarded as a general introduction to the subject. The final 
chapter deals with some of the auxiliary electronic circuits employed in conjunction with 
counters. This constitutes a very extensive subject in itself, so that, in the space at his 
disposal, the author has been able to give only a very brief outline of the types of circuits 
which can be employed. The book also includes a fairly comprehensive list of references 
to original papers on ion counters. Summing up, Professor Korff’s book can be recom- 
mended to those readers who want a good general introduction to the behaviour of ion | 
counters. J. L. MICHIELS. 


Theoretical Aerodynamics, by L. M. M1tne-THomson, M.A., F.R.S.E.. Pp. 
xx+363. (London: Macmillan & Co., 1948.) 40s. 


The purpose for which this book is written is to give a sound mathematical basis of 
theoretical aerodynamics on which applications to practical problems of aircraft design and 
flight must rest and to state clearly the underlying assumptions of the mathematical theory 
so that there is no uncertainty on what can legitimately be deduced from the assumptions 
as a first approximation. : 

The book has 19 chapters. Chapters I and II are introductory and concerned mainly 
with general features of fluid motion. Chapters III to XII and Chapter XIV deal in a 
commendable manner with well-known subjects of theoretical aerodynamics, namely 
two-dimensional motion, rectilinear vortices, the properties of two-dimensional airflow past 
a circular cylinder, Joukowski’s transformation theory of two-dimensional aerofoils, thin 
aerofoils, induced velocity, aerofoils of finite aspect ratio, lifting-line theory, lifting-surface 
theory and wind-tunnel corrections. The treatment of the mathematical theorems and of 
the illustrative applications chosen is made as simple as possible to help a reader whose 
mathematical knowledge does not extend much beyond the elements of the differential and 
integral calculus. The appeal which the selected applications make to an aerodynamicist 
familiar with the subject will depend on the extent to which his needs have been appreciated 
by the author : but the illustrations chosen give clues for the solution of related applications. 

Chapter XV on subsonic and supersonic flow serves as an introduction to a subject on. 
which a great deal of theoretical investigation has been made during recent years. Present 
knowledge on the mathematical theory of compressible flow is sufficient to provide subject 
matter for a complete book : it is not to be expected that it could be treated adequately in a 
single chapter. Chapter XIII and Chapters XVI-XVIII on propellers, simple flight 
problems, moments and aircraft stability respectively are limited in scope but suffice to give 
a student an insight into some theoretical problems related to the performance and behaviour 
of aircraft. Chapter XIX is well worth inclusion for the benefit of a reader who has little 
or no acquaintance with vector methods of presentation. 

Exercises for a student are given at the end of each chapter: lists of references to 
relevant papers are unfortunately not given. It is perhaps surprising that a chapter on the 
mathematical theory of boundary-layer flow is not included. 

The book fulfils the purpose for which it was written: it can be recommended as a 
good introduction to theoretical aerodynamics. A. FAGE. 


hy 


CONTENTS FOR SECTION B 


ia PAGE 
io} Editorial : ; : : ; ‘ : ; é ‘ : rage 
™) Dr. G. D. Wassermann and Dr. E. Wotr. On the Theory of Aplanatic Aspheric 
| Systems. , : : : : ; : . s : : eee 
m™ Dr. W. Exrenserc and Mr. R. E. Sipay. The Refractive Index in Electron Optics 
and the Principles of Dynamics : : ‘ : : : ' oe) 
™% Dr. H. H. Hopxins. The Disturbance near the Focus of Waves of Radially 
. Non-Uniform Amplitude . ; ‘ : : : : , : POD) 
| Dr. R. F. Barrow and Mr. E. F. Catpin. Some Spectroscopic Observations on 
; Pyrotechnic Flames ‘ : ‘ ; : ; : : ~ 832 
Mr. A. A. Jarre, Dr. J. D. Craccs and Mr. C. BaLakrisHnan. Some Experiments 
on Photo-icnization in Gases. : : : ‘ : Se AyIe Sey, 
Mr. G. C. Wixiiams, Dr. J. D. Cracces and Mr. W. Horwoop. The Excitation 
j and ‘Transport of Metal Vapour in Short Sparks in Air erwin ei , 49 
Dr. A. VAN ITTERBEEK and Dr. W. van Donincx. Measurements on the Velocity 
of Sound in Mixtures of Hydrogen, Helium, Oxygen, Nitrogen and Carbon 
Monoxide at Low Temperatures : E ; : 5 ; : Pers 
Dr. A. ScuattaMacn. Ultrasonic Dispersion in Organic Liquids ‘ ‘ eh!) 
Reviews of Books . : : : ‘ : ; : 5 : : eee? 
Contents for Section A . s : : : : : : : ; er! 
Abstracts for Section A’. : : : : f . ; : : 2 


ABSTRACTS FOR SECTION B 


On the Theory of Aplanatic Aspheric Systems, by G. D. WassERMANN and 
E. Wotr. 


- ABSTRACT. Methods are described for the design of two aspheric surfaces for any given 
centred system, so as to achieve exact aplanatism. The practical application of the methods 
is illustrated by the design of a reflecting microscope. 


The Refractive Index in Electron Optics and the Principles of Dynamics, by 
W. ExrensBerRG and R. E. Sipay. 


ABSTRACT. Inview of mis-statements made in the literature, the origin of the refractive 
index in electron optics is discussed in some detail, and the uniqueness of an expression 
_ previously given is demonstrated. On this basis, some general properties of electron 
optics are investigated. , 

| A relation between ray direction and wave normal is obtained. Whereas the refractive 
index is unique in terms of the magnetic vector potential A, this itself is arbitrary to some 
extent. It is shown that A must, for purposes of electron optics, be chosen so as to satisfy 
Stokes’ theorem and that, if it does, no observable etfects result from the arbitrariness of A. 
An expression for the optical path difference is given in terms of the magnetic flux enclosed. 
The results are applied to a number of questions, viz. the differential equations for tra- 
jectories, the focusing properties of an axially symmetric field and the interference pattern 
produced by two converging bundles of rays which enclose a magnetic flux. 


The Disturbance near the Focus of Waves of Radially Non- Uniform Amplitude, by 
H. H. Hopkins. 


ABSTRACT. Lommel’s problem, the disturbance near the focus of waves of uniform 
amplitude, is generalized to consider waves of radially non-uniform amplitude. Lommel’s 
numerical results are extended, and his U and V functions, tabulated on a new basis, 
together with some new related functions, are employed to provide extensive intensity 
distribution curves for the diffracting patterns associated with waves of both uniform and 


non-uniform amplitude. 
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Some Spectroscopic Observations on Pyrotechnic Flames, by R. F. Barrow and 
E. F. CaLpin. 


ABSTRACT. The results of a spectroscopic study of the white or coloured flames from 
typical pyrotechnic compositions are reported. Most of the visible light from these sources. 
comes from a limited number of atomic or diatomic emitters.’ The effects of secondary 
emitters on the dominant wavelength, colorimetric purity and relative intensity are con- 


sidered in general terms and some suggestions for improving the quality of the flames — 


are put forward. 


Some Experiments on Photo-ionization in Gases, by A. A. Jarre, J. D. Crass 
and C. BALAKRISHNAN. 

ABSTRACT. From an investigation of the discharge spread in Geiger counters filled 

with hydrogen, neon, argon or helium, it is possible to show that in certain circumstances. 

the discharge spread is due almost entirely to photo-ionization in the gaseous counter 

filling. The results enable the absorption coefficients of the operative radiations to be 

found. 


The Excitation and Transport of Metal Vapour in Short Sparks in Air, by 
G. C. Wituiams, J. D. Craccs and W. Hopwoon. 


ABSTRACT. The paper describes a study of the excitation temperature in certain spark — 


discharges of accurately known current characteristics. ‘The excitation temperatures are 
found by measuring intensity ratios for certain spectral lines where the relevant transition — 
probabilities are known, and are compatible with earlier work by Ornstein and his colla-_ 
borators on arc discharges. 

The excitation temperatures are, as would be expected, higher than those generally 
found for are conditions. 

Certain peculiarities relating to the evaporation of electrode metal are described. In 
particular the presence of discontinuous evaporation was noticed, and it was found that 
different electrode metals showed considerable variations in their behaviour in this respect. 


Measurements on the Velocity of Sound in Mixtures of Hydrogen, Helium, Oxygen,. 
Nitrogen and Carbon Monoxide at Low Temperatures, by A. VAN ITTERBEEK 
and W. van Donincx. 


ABSTRACT. Values of the velocity of sound in H,-N., H,-O., H,-CO, He—Ng, 
He—O, mixtures of varying proportions have been determined by experiment and 
compared with the theoretical values. Curves showing the variation of sound velocity 
with the pressure of the mixtures fo1 different temperatures (65° K. to 90° K.) are given. 


Ultrasonic Dispersion in Organic Liquids, by A. SCHALLAMACH. 


ABSTRACT. ‘The sound velocities v have been méasured, at frequencies f between 
about 700 and 3,000 kc/s., and at temperatures between about +15° and —80°c., of the 
following liquids: geraniol, di-dihydrocitronellyl ether, castor oil and iso-butyl alcohol. 
There is evidence of a positive dispersion (0v/@f>0) in the first three liquids at low 
temperatures. A negative dispersion has been observed in geraniol and di-dihydrocitron-- 
ellyl ether at higher temperatures, and in iso-butyl alcohol over the whole temperature 
range. 

The positive dispersions are discussed on the basis of the classical hydrodynamic theory, 
and the implications of a complex viscosity are examined. 
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